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1. Page 28, last paragraph. The second sentence should be changed to read "An application of the 
flux integral reveals further that Mq in Equation 94 becomes zero in this case, as is demon- 
strated in Appendix C." 

2. Page 29, Equation 96 should read 


n “1 


I< “ > ('•».) T F «)/ , . »», t. 

a. = 1 



1 ~ k T 

e 


N 

“>l ( +k a) P / (Mi) 


r 0 


-r/M 0 


1 + Mi/Mo [j 


N 

- 





S, ^ Si t 

k) p - M 

+ 3 M 0 M n r ( S 0 “ 1 

; i = ±1, 

•, ±n) . 

(96) 


The foregoing changes have become necessary because of an incorrect derivation in Appendix C. 
The following is a corrected version of the relevant derivation. 


3. Page 70 et seq. The section from Equation C9 through Equation C21 should read 
Upon substituting A. = 1 in Equations C8 and C9, Equation C6 reduces to 
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Upon setting X. = 0 in Equation C7, and remembering that ~ 0 > Equation ClOd further reduces 
to 


= fJ Mo e" TAo + (l “ 25 0 ) 


F(t) = F 0 <^ 0 




+ k„T “k_T\ -vVo 

e - M a e j - /j - 0 T 0 e 


l_ a«l 


(Cll) 


The net flux crossing the top of the cloud is found, upon setting r = 0, from 
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and crossing the bottom of the cloud, upon setting r = r t , from 
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If the cloud is semi-infinite, all M_ a = 0 ( a = 1 , • • - , n ), and Equation C12 becomes 
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In the special case of conservative scattering (£ 0 = l), Equations C12-C14 are no longer 
necessarily valid relations. In this case the azimuth- independent term of the intensity in the case 
of a finitely thick cloud is given by (Equation 94) 
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where M ±a ( a = 1, ■ * • , n - 1), M 0 , and M n are the 2 n constants of integration. It is clear from 
Equation Cll and the preceding discussion that in this case the flux integral must reduce to 
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Equation C16 reduces after some algebra to 
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This is the net flux of the diffuse radiation field. Since the net flux of the direct field (from 
the point source) is given by 
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the net flux of the total radiation field is given by 
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The flux integral F T (r) is seen to be independent of optical depth, as it must be for conservative 
scattering. 

If the cloud is semi-infinite and 2> 0 = l, then, as r - 0, Equation C18a must reduce to 

F(0) = m„f 0 : (C19) 


i.e., the "albedo" of the cloud must be unity. Thus M 0 is zero, and the expression for i<°> (t,/a.) 
in this case is seen by inspection to be 


lC°>(T, Mi ) = 


4 F o 



1 -k _ r 

M a e ° 

1 + Mi k a 


N 

€l ( +k a) p ; (Mi) 


“ r Ao 

r 0 e 

+ 1 + Mj/Mo 


N 

Z] ft (-ik) ft N 


+ 3(i Q M n 


\-i - o 


(C20) 


where M a (a = l , ■ • • , n - l) and M n are the n constants of integration to be determined from 
the n boundary conditions 


I<°> (0,-Mi) = 0 ( i = 1, n) . 

Continue with the section on The Law of Darkening, page 72. 


(C21) 


4 


Issued 9-2*f-65 


NASA-Langley, 1965 



TECH LIBRARY KAFB, NM 



)Db? c JbS 

NASA TR R-215 


RADIATIVE TRANSFER IN A CLOUDY ATMOSPHERE 

By R. E. Samuelson 

Goddard Space Flight Center 
Greenbelt, Md. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Office of Technical Services, Department of Commerce, 
Washington, D.C. 20230 -- Price $3.00 



RADIATIVE TRANSFER IN A CLOUDY ATMOSPHERE 


by 

R. E. Samuelson 

Goddard Space Flight Center 

SUMMARY 

The equation of radiative transfer in the context of a partially 
thermally emitting, partially anisotropically scattering plane -parallel 
cloudy atmosphere is derived. The derivation allows an exact inter- 
pretation of the auxiliary quantities in the equation of transfer in terms 
of the Mie scattering parameters. Explicit solutions are given by the 
method of discrete ordinates in accordance with Chandrasekhar’s pro- 
cedure, and extended to include thermal emission at infrared wave- 
lengths. Solutions of this type, restricted to plane -parallel layers 
bounded on both sides by a vacuum and characteristic of a phase function 
for single scattering independent of optical depth, are referred to as 
solutions to the restricted problem. To extend the treatment to the 
general problem, a procedure wherein layers with different scattering 
properties are combined is developed in terms of the restricted solu- 
tions; two explicit examples are worked out. Such quantities as the 
angular distribution of outgoing radiation and the net outgoing flux are 
consequences of these solutions. 
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RADIATIVE TRANSFER IN A CLOUDY ATMOSPHERE 


by 

R. E. Samuelson 
Goddard Space Flight Center 


INTRODUCTION 

The study of radiative transfer in a cloudy atmosphere is extremely involved due to the fact 
that clouds are composed primarily of particles having physical dimensions comparable to the 
wavelengths of radiation of interest. This realm of particle sizes requires a much more elaborate 
treatment of the theory of single scattering (scattering of radiation by one particle) than do the 
theories for single scattering for particles of quite different sizes. 

Even where the theory for scattering by single particles is well understood, it still remains 
to describe the process of the transfer of radiation through a field of these particles, and, in par- 
ticular, to describe the angular distribution of radiation from such a field. It is due to the order of 
difficulty in describing this compound process quantitatively that many investigators in the past 
have been led to treat various simplifications of the theory. An example of one such simplification 
involves adding a diffuse (isotropic) component to the intensities derived from a regular single 
scattering problem. Another is the two stream method used by various authors, whereby the es- 
sential assumption made is that the radiation field may be divided into two streams normal to a 
stratified cloud layer. Havard (Reference l)has improved upon this method by replacing the specific 
intensity with the flux density as the dependent variable. Both methods suffer from the difficulty 
that nothing can be said about the angular distribution of radiation emitted from the cloud. 

On the other hand, Churchill etal. (Reference 2) have obtained exact solutions for the angular dis- 
tribution of outgoing diffusely reflected and transmitted radiation which are valid for three-term 
single scattering phase functions, and have developed a program which will solve transfer problems 
involving phase functions of much greater complexity exactly, provided among other assumptions 
that the phase function is independent of optical depth. The major limiting factors are the number 
of terms carried and the time involved on the electronic computer. A further limitation, which in 
practice is relevant only at wavelengths of a few microns and beyond, is that no account of thermal 
emission can be made. 

It will be the purpose of this paper to formulate in a logical manner the appropriate equation 
of radiative transfer consistent with restricted physical situations, indicate how numerical solutions 
may be obtained for these situations, and finally, show how these restricted solutions lead in a 
natural way to more general solutions for physically realistic model cloudy atmospheres. In order 
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to interpret the various single scattering parameters contained in the equation of transfer, it is 
convenient to derive this equation in the context of some single scattering theory. We have chosen 
to adopt the Mie theory implicitly for this purpose. 

Beginning with the fundamental concept of individual photon-particle interactions, we are led 
in a natural way to the exact interpretation of such physically significant quantities as the phase 
function for single scattering, the albedo for single scattering, and the optical depth, all of which 
are explicitly contained in the basic equation of transfer. This development further illuminates 
the nature of the various approximations which are employed in order to make the scalar equation 
of transfer, in the context of polydispersed plane -parallel media, amenable to solution in some 
finite order of approximation. 

Once the equation of transfer is formulated we proceed to the results of an analysis by 
Chandrasekhar. In particular, the solutions for the angular distribution of diffusely reflected and 
transmitted radiation at the surfaces of a plane-parallel cloud are given. In addition, an extension 
of this analysis to include radiation thermally emitted by a cloud having an arbitrary (but known) 
temperature profile is made, and analogous solutions to the more extended problem are derived. 
The solutions in all cases are obtained by the method of discrete ordinates compatible with the 
relevant equation of transfer. These solutions are referred to as solutions to the "restricted" 
problems; i.e., the relevant strata are assumed to be bounded on both sides by a vacuum, and are 
irradiated at most by a single outside point source. 

The problem is then extended to include the effects of the surrounding atmosphere and ground. 
The two explicit problems considered involve the solutions for the net flux and angular distribution 
of: 1) radiation which is diffusely reflected from the top of an optically thin atmosphere overlying 
an optically thick cloud, and 2) radiation which is thermally emitted from the top of an atmosphere 
containing an optically thin cloud. In the latter case only the cloud is assumed to contribute to the 
scattering of radiation. 

The formulation of the first of these more general problems consists of a linear integral 
equation containing as free parameters solutions to the more restricted problems previously dis- 
cussed. This integral equation is shown to reduce to a system of independent linear integral equa- 
tions of the second kind containing one independent variable only. Once these equations are solved 
by standard methods the remainder of the solution becomes quite straightforward. 

The formulation of the second of these general problems consists only of a combination of the 
solutions, and integrals of these solutions, of the relevant restricted problems previously alluded 
to. By virtue of formulating the problem in the context of only one scattering layer, no integrals 
containing the dependent variable appear in the formulation. The other layers (including the ground) 
are assumed to contribute only purely emitted radiation. 

Many extensions to the basic theory are possible. The solutions obtained here are primarily 
restricted to those which are not excessively difficult to solve numerically with the aid of an 
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electronic computer. Inherent in the final solutions are any defects in the formalism due to the 
following nine major assumptions: 

1. plane-parallel atmospheres, 

2. unpolarized radiation, 

3. spherical cloud particles, 

4. well- separated cloud particles, 

5. phase function for single scattering independent of optical depth, 

6. local thermodynamic equilibrium, 

7. blackbody emission from the ground at long wavelengths, 

8. finite Legendre polynomial expansion of the phase function for single scattering, 

9. finite approximation solution to the equation of transfer. 

Besides these approximations, there is the additional restriction of having to approximate functions, 
which vary continuously with certain physical parameters, by step functions. Examples are the 
variation of the complex index of refraction with wavelength and the variation of particle distribu- 
tion with particle size. It is true that these steps may be taken rather fine, but practical consider- 
ations of time and efficiency are going to limit the accuracy. 

In the above approximations, it is felt that only the concept of plane -parallel atmospheres and 
the independence of the phase function for single scattering upon optical depth are subject to severe 
criticism. A partial compensation can be made for polarization in the sense that those intensities 
which are incorporated in the general solutions which result from a Rayleigh scattering atmosphere 
will be considered to be those intensities which result from the correct matrix formulation of the 
equation of transfer for the relevant restricted problem. The other approximations are not felt to 
be particularly severe. Perhaps the concept of spherical cloud particles should be modified; cirrus 
clouds are an example of one physical situation which might require this. If the cloud were divided 
into more than one layer, and each layer treated separately, the approximation of requiring the 
phase function for single scattering to the independent of optical depth could also be modified. 
However, this would make the general problem ,that of combining the separate restricted solutions 
into a general solution through the use of certain integral relations to be derived in the last section 
of this paper, extremely complex, and hence will not be considered. 

Five appendices are included in this paper. 

Appendix A shows that the phase function for single scattering is normalized to the albedo 
for single scattering when this phase function is approximated by a finite series expansion of 
Legendre polynomials. 

Appendix B gives the solution for a particular integral satisfying the equation of transfer which 
includes effects of thermal emission. 
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Appendix C gives solutions for the various flux integrals which arise from solutions to the 
equation of transfer. 

Appendix D suggests methods of numerically solving certain integrals and integral equations. 
The integrals are replaced by Gaussian sums, and each integral equation is correspondingly re- 
placed by a system of n equations in n unknowns. 

Appendix E contains a list of symbols used in this paper. 


FUNDAMENTAL QUANTITIES 
Introduction 

In a gross way we may think of observables, namely the radiative specific intensity, outgoing 
flux density, and so on, as being field quantities, and we will want to inquire how clouds affect the 
radiation field as a whole. But in order to do this it will be necessary to investigate individual 
photon- particle interactions, for these individual processes are responsible for variations of the 
field. 

Regardless of the physical processes involved, the study must be undertaken in some frame of 
reference. The kinds of clouds which will be adopted as models will be plane- parallel, i.e., strat- 
ified with no horizontal inhomogeneities. It can be anticipated that the optical path length of radia- 
tion in clouds is very short compared to the radius of curvature of any planetary atmosphere, and 
that vertical inhomogeneities will far outweigh those of horizontal extent. Only ’’corrugated tops” 
of clouds are expected to cause serious difficulty; in fact this deviation from plane-parallelism is 
in practice probably the most serious defect of the models. 


The Coordinate System 

Since it will be necessary to consider such processes as absorption, scattering, and emission 
of radiant energy on a microscopic scale, we shall need the concept of an element of mass dm con- 
taining many particles, restricted to a volume dV, and characterized by a mean density p of matter 
in the element. The position of this mass element will be specified by the vertical distance z meas- 
ured positively from the ground to dm. Directions at dm will be specified by the cosine of the zenith 
angle e> denoted by p, and the azimuthal angle 0 (Figure 1). The zenith angle e is measured posi- 
tively from zero (the zenith) to v (the nadir). The azimuthal angle 0 is measured through 2 n radi- 
ans in the plane of stratification from some arbitrary angle 0 O . Both the values of 0 O and the 
direction of measurement will be specified at a later time. In general directions will be indicated 
by the symbol (p , 0). For example, i(z, /x, 0) is the intensity of radiation at dm (at a level z) in the 
direction (0, 0). This direction is inclined to the normal to the plane of stratification by an angle 
6, and is contained in the plane defined by the normal and the azimuthal angle 0 . 


4 



The Specific Intensity 


The amount of radiant energy which is 
transmitted across an elemental surface area 
dcr in a time dt and in a direction inclined at an 
angle 0 to the normal to dcr will be given by . 
SE v is further restricted to the frequency range 
(v, v + d^) and to a solid angle dco. The depend- 
ent variable of principal interest is the "mono- 
chromatic” specific intensity defined by the 
limiting ratio 

cos 0 dcr dcodv dt (i) 

where da, do >, dv , dt -* 0 in any manner. Thus i v 
is the rate at which radiant energy confined to a 
unit solid angle and unit frequency interval 
crosses a unit surface area which is normal to 
the direction of radiation. 

The Phase Function for Single Scattering 


Z 



Figure 1 — The relation between the scattering angle 0 
and the coordinate ang les <p, <p‘ , 0, and 0 '. Radiation is 
assumed to be incident on the mass element dm at 0 in the 
direction 0P 1 = ( f.i , <p) and scattered by dm through the 
angle© into the direction 0P 2 


We shall need the concept of a phase func- 
tion for single scattering, fP(cos®), which de- 
scribes the angular distribution of radiation 

scattered once. If AE v represents the fractional amount of energy E v incident on dm in unit time in 
the direction (^, <p), which is either absorbed or scattered in all directions in unit time, and d(AE„) 
is that fraction of AE v which is scattered into the direction (/x', $') and contained in the solid angle 
dco ' , we may formally represent this fraction of scattered radiation by 


d[AE„ ( z - m' . 4 >‘ )] 


5>(cos®) AE„ (z, 


dco ' 

M. <t > ) 4tt 


(2) 


where © is the angle through which the radiation is scattered.* 


If the albedo for single scattering is defined to be the ratio of radiation scattered in all di- 
rections in unit time to the radiation extinguished (absorbed plus scattered) in unit time, we clearly 
have 


AE V (z, /_i, 


J d[AE„ (z, fj.' , </>' )] - co 0 - J , :P(cos ®) ^ 


( 3 ) 


*The. tacit assumption has been made that multiple scattering does not take place in dm. The consequences of this assumption will be 
made apparent in the next section. 
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where the integration is performed over all solid angles. The function ?(cos ©) may formally be 
replaced by 


fP(cos0) - pO',0'; /x, 0) (4) 

since obviously the scattering angle can be expressed in terms of the coordinate system; 
p(/x', 0' ; m, 0) then refers to radiation originally in the direction (/x, <p) which has been scattered 
into the direction , 0'). If 2> 0 = i we have conservative scattering (no absorption), and 
p(/x f , <p‘ ; fx , 0) is normalized to unity. 

The relation between cos 0 and 0? mS and 0' can be easily obtained from Figure 1. Let us 
suppose that radiation is incident on a mass element dm at 0 in the direction 0P p and is scattered 
at 0 through an angle © into the direction 0P 2 . From the spherical triangle ZP 1 P 2 we obtain the 
following relation: 


cos© “ cos 0 cos 6' + sin 6 sin 6 ' cos (<fc‘ - 0) - jx jx' + (l - /x 2 ) / (l~^' 2 ) / cos(0'-0) (5) 

where fx - cos 6 and ^ - cos 6‘. The angles 6 and O' are respectively the zenith angles of incidence 
and scattering, and 0 and 0 ' are the corresponding azimuthal angles. It should further be noted 
from Equation 5 that p(/x', <p f ; jx, <p) obeys Helmholtz ! s principle of reciprocity; i.e., 


p(m'»0'; m. 0) " p(m» 0; 0',0') ■ 


(6) 


The Net Flux 

Equation 1 gives the rate at which energy confined to a unit solid angle and unit frequency in- 
terval flows across a unit area which is normal to the direction of radiation. The "monochromatic" 
net flux tt¥ v is the net rate at which energy per unit frequency interval flows across a unit surface 
in all directions, and this is given by 


= 

where the integration is to be performed over all solid angles. It should be noted that Equation 7 
gives, in effect, the difference between the upward and downward fluxes across a unit surface area; 
therefore, deviations from a null net flux reflect deviations from conditions equivalent to those in- 
side a perfectly insulated isothermal cavity. 

In terms of the coordinate system previously specified the net flux becomes 




T cos 6 doj 


(V 


77 F 


V 





fx, <fi) d/xd<fi . 


( 8 ) 
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It is clear from Equation 8 that ¥ v is in general a function of z . An extended account of the explicit 
integrations which will be required in the context of the specific problems encountered in this 
paper is given in Appendix C. 


BASIC CONCEPTS AND THE EQUATION OF RADIATIVE TRANSFER 
Introduction 

It will be easiest to consider the interaction of radiation with matter from the "Lagrangian" 
point of view; i.e., the movement of individual particles (photons in our present context) will be 
followed. However, the observable parameter of interest, namely the monochromatic specific in- 
tensity i v , is a field of these photons, and we shall be interested in following variations of this 
field as it interacts with matter in some specified way. For this purpose it will be more practical 
to study local variations of i v without regard to the individual history of each photon. This is the 
"Eulerian" point of view (Reference 4). 

An abrupt transition from the Lagrangian to the Eulerian point of view is somewhat unappeal- 
ing in that it seems rather artificial. In order to appease this feeling of discontinuity of logical 
consequence we shall endeavor to smooth out the transition by retaining something of the Lagran- 
gian point of view while introducing the Eulerian point of view, and in doing so give some kind of 
picture illustrating the essential identity of these two concepts. 

Classification of Photon-Particle Interactions 

Consider the photons of frequency v interacting with dm to be classified by "type", each type 
depending upon the result of interaction as well as upon the intrinsic characteristics of the photon 
itself. Each class of photons is therefore a field of these photons which behaves exactly as every 
member of the field individually. We will want to restrict our attention to one field at a time, which 
in essence is the same as restricting our attention to one photon of this field at a time. Thus, the 
field from the Lagrangian point of view (in a less strict sense of the phrase) will be followed. 

Since each photon behaves (by definition) like every other photon in the field, it is clear that 
not only must each photon be identical with every other photon of this class, but also the system of 
particles which this class as a whole interacts with must be composed of exactly the same kind 
(dimensions, refractive index, etc.) of individual particles in order that the separate interactions 
be identical; and if more than one interaction per photon takes place, the order, number, and 
character of these interactions must be the same for each photon. Since the last restriction 
creates an insurmountable strain on the imagination, it is seen immediately that we must restrict 
the volume dv of the mass element dm to dimensions considerably less than the mean free path of 
an individual photon, so that only single interactions are possible in dm. The additional problem 
of correctly specifying the orientation of each particle can be circumvented only by limiting our- 
selves to homogeneous spherical particles. 
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We define four classes of photons, each photon in the frequency (energy) range {v, v + dv), and i 

each interacting with one (spherical) particle in the radius range (r. , r i + dr i ) which is of homo- 
geneous composition of complex refractive index n - i*. The four classes are distinguished by 
four different types of interactions, and are defined as follows (cf. Figure 2)*: | 

1. That system of photons which is incident upon dm in a time dt and in the direction (u, 0 ) 
contained in the solid angle do>, and is (singly) scattered into the solid angle da/ in the direction 

4>') by interactions in dm with particles in the radius range (r £ , r A + dr.). This scattering 
process may either be considered as a scattering of a certain fraction of the number of incident 
photons into da/, or as a re-direction of a fraction d[SE v (z, fi‘ , <f >‘ )] of the incident energy 8 e v (z, /z, 0) j 

into do/ . 

2. That system of photons which is incident upon dm in time dt and in the direction (^, 0 ) con- 
tained in the solid angle d^, and is absorbed by the particles under consideration in dm. This ab- 
sorption process may be considered as an absorption of a certain fraction of incident photons or 

as a diminution of the incident energy 8e v (z, /z, <p) by an amount d[SE v (z, /z, 0)] . I 

3. That system of photons which is incident on dm in time dt and in the direction (^', 0') con- 
tained in the solid angle da/, and is (singly) scattered into the solid angle doj in the direction (j±, 0). 

Again, this scattering process may be considered either as a scattering of certain fraction of in- 
cident individual photons into dco 7 or as a re-direction of a fraction d[&E v (z, /lz, 0)] of the incident ; 

energy 8E V (z, p , 0') into doz 

4. That system of photons which is (thermally) emitted from the particles under considera- 
tion in dm into the solid angle doj and in the direction (u, 0 ) in a time dt . This thermal emission 
process may be considered either as an emission of many individual photons, each in the energy 
range (hv, hv + dhv ), or as a source of energy d[SE v (z, /z, 0 )] which is emitted into do; in the direc- 
tion (/z, 0). In principle we need not have considered only thermal emission; however, in planetary 
atmospheres we would not expect other types of emission to play a very significant role. 

The first two classes of photons are those which are lost from the radiation field in the direc- ; 

tion (tz, 0 ) by scattering and absorption respectively. The last two classes consist of photons which \ 

are gained by the radiation field in the direction (/x, 0 ) respectively by scattering and emission. j 

These are obviously not the total losses and gains to the radiation field, however, since only inter- 
actions with particles in the radius range (r . , Ti +dr i ) have been considered. Later on an integra- . \ 

tion over all particle sizes will have to be made. *; 

The four classes of photons have the following characteristics in common: j 

a. frequency range (v, v + dv). j 

b. interaction with one (spherical) particle in the radius range (r. , r. + dr. ). 

c. interaction with particles of homogeneous composition of complex refractive index (n - i*). i 

\ 

1 

j 

*Note that do> and dco are elements of solid angle which are both subtended at dm. j 
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The classes are also characterized by the differences tabulated below: 


Class Characteristic 

Incident 

Direction 

Solid Angle 
Subtended at dm 

Emergent 

Direction 

Solid Angle 
Subtended at dm 

Singly Scattered (Loss) 

(m. <t>) 

d co 

(m’> <P') 

do/ 

Absorbed (Loss) 

(m, <P) 

d co 

None 

None 

Singly Scattered (Gain) 

(m',** ) 

dto' 

(m, 4>) 

dco 

Thermally Emitted (Gain) 

None 

None 

(m, <£) 

d co 


The classes may be characterized in terms of energy as follows: 

1. Re-direction of fraction d[§E v (z, ^ , <£')] of incident energy 8E V (z, /x, <p) into d V. 

2. Absorption of fraction d[SE^ (z, /x, 4>)] of incident energy SE V (z, /x, <f>) by dm. 

3. Re-direction of fraction d [se v (z, /x, <£)] of incident energy SE V (z, /x' , <£' ) into dco. 

4. Emission of energy d[§E v (z, /x, <£)] into do>. 

Single Scattering 

In order to describe the interaction of one photon with one particle quantitatively, an appeal 
will have to be made to some single scattering theory. The Mie Theory (Reference 1) will suit our 
purpose if we restrict ourselves to homogeneous spherical particles in a radiation field of plane 
waves.* Conceptually the Mie theory describes the outgoing "scattered" radiation in terms of a 
spherical wave front which depends upon the wavelength of the incoming plane wave and the size 
and complex index of refraction of the spherical particle. The radial component of the outgoing 
spherical wave dies out as the square of the distance, while the tangential component dies out as 
the first power of the distance. Thus, the outgoing spherical wave tends to take on more and more 
the character of a plane wave as the radial distance from the scattering center is steadily in- 
creased. Since it has been implicitly assumed that many of the incident plane wave "photons" have 
originated from scattering and emission processes in other mass elements, it is clear that the di- 
mensions of the mass elements must be sufficiently large such that the outgoing spherical wave of 
scattered radiation from a particle in element dm', say, has traveled a sufficient distance to be- 
come essentially a plane wave by the time it suffers another interaction with a particle in dm. In 
a crude way this means that the mean free path of a photon from mass element to adjacent mass 
element must be large compared to the wavelength of radiation, which is going to place an upper 
limit on the size of the particles compared with their mean separation distance. We shall always 
assume that physical conditions prevail such that the Mie theory gives a valid picture of single 
scattering in clouds. 

It will not be the purpose of this paper to describe the Mie theory. Van de Hulst (Reference 
6) and Born and Wolf (Reference 5) give a good account of it, and Havard (Reference 1) gives rather 


♦Photons are now regarded as plane waves, with the realization that this rather abrupt transition of our "picture* of a photon leaves 
something to be desired. 
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complete details on how to compute the phase function for single scattering [?(cos©) = pO, <f>; ^ , <£')] 
and the efficiency factors* for extinction, absorption, Q A (i) , and scattering, Q s (i) . The relation 

Q E (i) = Q s (i) + Q A (i) (9) 

always holds. With this brief digression we return to a more detailed discussion of the first class 
of photons. 


The Radiation Field; Losses and Gains 

Construct a convex closed surface SS around dm = pdv such that the volume contained by SS is 
large compared with dv but small otherwise (Figure 2). Let dA and dA' be elements of SS such that 
the direction from d A to dm is (^, <t>) and the direction from dm to dA' is (//, <£'). Further let d co and 
d a*' be respectively the elements of solid angle containing dA and dA' as seen from dm, and da and 



Figure 2 — Illustration of the elements of solid angle da; and dad subtended at the mass element 
dm respectively by dA and dA ' , elements of the convex bounding surface SS. Also shown are 
the elements of solid angle da and da', subtended respectively at dA and dA ' by dm. 

♦Each efficiency factor is defined as the ratio between the effective cross-section of the particle to the geometric cross-section 
7tr. 2 , where r. is the radius of the i th particle. Thus the efficiency factors for extinction, absorption and scattering are 


q e ( ° 

q a <‘> 


Q (i) 

v s 



( X ^TTX 


0 
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da' the elements of solid angle containing dm as seen respectively from dA and dA'. Referring to 
Equation 1, it is seen that the amount of energy SE V (z, /z, <£) crossing dm in a time dt and in the 
frequency range v + dv ), which has originated outside SS and has also crossed dA, is 


SE V (z, (jl, <£) - l v (z , /z, 4>) Mi dA da dv dt 


(10) 


where fj. 1 is the cosine of the angle between the direction (/x, 4 ) and the normal to dA. 

Of all the energy crossing dA which is contained in the solid angle da, a certain fraction will 
be singly scattered by dm into the solid angle do/. The fraction of energy which is thus scattered 
in a time dt is seen to be (Equation 2) 


d[8E„ (z, m', <£')] 
SE V (z, /x, <t > ) 


P (i) P 

r S *1 


da/ 

( m ' . 4 >' ; <£) 477 


( 11 ) 


where P s (i) is the probability that any one photon which is incident on dm is scattered in any direc- 
tion whatsoever; p i Ou.', </>'; ^ <f>) in this context is the phase function for single scattering through 
the angle between the directions Oz, <$>) and (^', <£') normalized to unity. Now the probability of 
scattering p <*> of one photon is just the ratio of the total available effective scattering cross- 
section of all the particles in the radius range (r i( r . + dr i ) contained in dm to the geometrical 
cross-section da of dm as seen in the direction (^, <f>), and this is given by 


P (i) 

s 


Qs (i) 77r. 2 Nj dr. dV 


da 


(12) 


where dv is the volume of dm, n* is the number of particles per unit volume per unit radius range 
centered about r . , and Q s (i > is the efficiency factor for scattering from the Mie theory as described 
previously. Equation 12 is valid so long as there is no "shadow” effect, i.e., the probability that 
any one particle is screened off from any other particle is negligible; this requires that only single 
scattering prevails, or, put another way, that P s (1 > is small compared to unity. Collecting Equa- 
tions 10, 11, and 12 we have: 


d [SE v (z, fi 1 , <j> 9 )] 


Q <*> 77 r . 2 N. dr x dV 

'do- 1 “ “ l v (Z, M. 4>) Ml dA da di^ dt p. (/z # , <j>' ; 


d co 1 

M. <£) 4^ 


(13) 


There must clearly be a loss of intensity d[s_ i s (z, /z, <£)] in the direction (jj. , 4>) which is as- 
sociated with the loss of energy d[8E„ (z, m' , <t >' )] from the direction , <f>), and this is given by 
(cf. Equation 10) 


d [SE V ( z , /x' , <£' )] - d [S_ I s (z, /x, <£)] /ij dAda dv dt 


(14) 
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Comparing Equations 13 and 14, we obtain 


d[s_ I s (z, /i. <£)] 


Q s « 7rr . 2 Nj dr. dV 

^ I v (z, /x, <£) Pi ( M ', d>'; 


doj ' 


(15) 


Equation 15 is valid for scattering of incident photons in the frequency range (v, v + dv) into 
the solid angle dco' by particles in the radius range (r. , r. + dr. ) . In order to obtain the effect of 
scattering from particles of all sizes in all directions, Equation 15 must be integrated over all r. 
and a/. Thus, the total intensity in the frequency range (v, v + dv) and in the direction 0) which 
is lost from this direction by scattering from dm is 


IN 0 dV du) > 

S_I s (z,M, <t>) = l v ( z > 0) Q s (i) 7TT i D i Pi O' , 0' ; 4>) dr. , (16) 

Jo Joi‘ 

where N 0 is the total number of particles per unit volume and (the relative number of particles 
per unit radius range centered about r A ) is the normalized distribution function of particle sizes; 
i.e., N. - N oDi . 

By definition ( i/4tt) p.(/i', 0 '; 0 ) integrated over all solid angles is unity regardless of the 

particle size involved; however we shall delay this inter gration until later. 

We turn our attention now to the second class of photons defined previously. From Equations 
10 and 11 and the corresponding discussion it should be clear that the energy in the frequency 
range (v, v + dv) crossing dA (Figure 2) in a time dt, which is absorbed by the particles in the 
radius range (r., r. +dr i ),is given by 

d [§E V (z, fjL, 0)] - l v (z, fjL, 0) fx l dA da dv dt P A (i > , (17) 

where P A ci >, the probability of absorption of a single photon in the energy range (hv, hv + dhv), is 
given by (cf. Equation 12) 


P« = 


Q a <‘> vt? Nj dr. dV 


do- 


(18) 


By analogy with Equation 14, the loss d [&_ i A (z, fi, <#>)] of the intensity I v (z, i< , i ) which is incident 
on dm in the direction (jx, d>) is given by 


d [SE„ (z, fi, *)] - d [8_ I A (z, fi, <£)J jLtj dAda dv dt 


(19) 


Comparing Equations 19 and 17, we obtain, with the aid of Equation 18, 


Q a <‘> rrr^Nj dr t dV 

d[S_ I A (z, n, </>)] = 1 l v (z,n,<j>). 


( 20 ) 
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and this integrated over all particle sizes becomes 


S_ I A (z, /X, 0 ) 


n 0 av r 

— 35- I v (z, /x, 0) Q A (i) 7TF. 2 D a dr. . 
^0 


( 21 ) 


Equation 20 is the total intensity in the frequency range (v, v + dv) and in the direction (^, 0) which 
is lost from this direction by absorption in dm. 

In order to establish the gain to the radiation field by the scattering of the third class of pho- 
tons into the direction (^, 0), it is necessary only to reverse the sense of direction of the radiation 
field in Figure 2. By once again tracing out the consequences of the scattering process it becomes 
clear that Equations 10 to 13 remain valid upon an interchange of the primed and unprimed quan- 
tities. Thus, the gain of energy d[SE v (z, 0)] in a time dt and in the direction (a, 0 ) contained in 

the solid angle dco, which has resulted from a scattering of the energy by particles in the radius 
range (r. , r i + dr.) in dm, and which was originally in the direction (/x', 0') and contained in the 
solid angle da/, is 


r Q s (l) 7Tr i N- dr { dV 

d[SE v (z, /x, 0 )] = foT> ' l v (z, /x' , 0 ' ) /x/ dA' da' dv dt p £ (/x, 0 ; , 0 ' ) * 


d " (22) 


Here da-' is the geometrical cross-section of the mass element dm as seen in the direction (/x', 0'), 
da' is the solid angle subtended by do-' at dA', and /xi ' is the cosine of the angle between the direc- 
tion (fj,' , 0 ') and the normal to dA' . 

The corresponding gain in intensity at dm in the direction (^ , 0) must be given by 


d [SE V (z f fj,, 0)] - d[S + I s (z, /x, 0)] do-doodvdx. , 


(23) 


since d[SE v (z, f±, 0)] is just the energy in a time dt which crosses normal to the surface element 
do- at dm and is contained in the solid angle dco. Comparing Equations 23 and 22, we obtain 


Q s (i) rrr 2 N. di^ dV ^ dA' da' 

d[S+ I s (z, /x, 0)] = 3? I v (z, /x\ 0') p i (/x, 0; fj.' , 0') 4?Tdcr ' 

If the distance between dm and dA' is denoted by r, then, from the geometry, 

fjL t ' dA' = T 2 dco' 

and 

do - ' = r 2 da' . 

With the aid of Equations 25 and 26, Equation 24 becomes 

Qs^Ttr^N. dr. dV da/ 

d [S + ! s (z, 0)] = do- Pi (M. 0; /*'. 0') I v ( z > M'f <#>') 4 F * 


(24) 


(25) 


(26) 


(27) 
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In order to obtain the total contribution to the intensity ir the direction ( M , 0 ) by scattering from 
dm, Equation 27 must be integrated over all solid angles and all particle sizes r.; i.e. 

(cf. Equation 16), 


i>Q uv j | 

B + i s (z, M. 4 > ) = ~d5“ J Q s (i)77T i 2 Pi Pi (m. 0; M', <£') I v (Z, M # . <£') -35T dr. . (28) 

The fourth class of photons describes the contribution to the radiation field by thermal emis- 
sion from dm. In order to calculate this contribution, consider the surface of Figure 2 to be a 
perfectly insulated isothermal enclosure in which dm is maintained at a constant equilibrium tem- 
perature T. Consider N i dr A particles in the radius range (r., r. +dr A ) contained in dm, each of in- 
ternal temperature T. Since the radiation field within the enclosure is in equilibrium with its sur- 
roundings and isotropic, the amount of energy in the frequency interval (v, v + dv) and in the 
direction (^, <f>) contained in the solid angle dco, which would be emitted from dm in a time dt upon 
an instantaneous removal of the cavity walls, is given by (cf. Equation 23) 

d [8E V (z, /X, 0)] = (T) dr A x A (i) dVda;dvdt , (29) 

and this must be equivalent to 

d [§E V (z, /X, 0)] = d [S + I E ( z , At, 0)J dcr dco dv dt , (30) 

where d[s + I E (z, {i, 0)] refers to the gain in intensity in the direction (ji , 0) due to thermal emission 
from dm, and where b v (t) is the Planck function in intensity units. If the total number of particles 
of all sizes per unit volume is n 0 , and the size distribution in dm is given as before by D i? then 
Equation 29 may be re-written as 


d [SE V (z, /z, 0)] - N 0 dVQ A (l) 77r. 2 Di B v (T)da;dvdt dr £ - N 0 dV B v (T) dw dv dt Q A (i) 77r. 2 D. dr. , (3f) 

J SE l/ Jo Jo 

where x A (i) = Q A i>77r i 2 * Relating Equation 30 to Equation 31 we have 

f N 0 dV 

S + I E (z,/i, <f>) - d[S + I (z, 4>)] = ~ cb - (T) 

If the enclosure is not replaced, dm will be subjected to the anisotropic radiation field of arbi- 
trary energy density. What happens now is largely a function of the relative importance of: (1) 
collisions between molecules contained in the particles, compared with (2) interaction between 
these molecules and the radiation field, as a cause of molecular absorptions and emissions. If 
interactions with the radiation field predominates, the emission will essentially consist of the 
spontaneous emission of photons from excited molecules, and induced emission of photons from 


j 


Q d) 77r 2 D. dr . 


(32) 
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excited molecules through perturbations due to the radiation field; the latter type of emission is 
proportional to the incident intensity and is therefore anisotropic in the same sense as the radiation 
field. The molecules will be excited through the absorption of incident radiation, and the local tem- 
perature will be strictly dependent upon the photon density. Thus the radiation emitted from the 
particles cannot be isotropic because of the contribution from induced emission, unless the radiation 
field itself is strictly isotropic. 

At the other extreme, collisions! battering of molecules predominates, and the thermal (iso- 
tropic) emission far outweighs emission induced by the radiation field. This will occur where the 
density of molecules is great enough so that the frequency of encounters among molecules is much 
larger than the frequency of encounters between molecules and incident photons. We shall suppose 
this latter condition always to prevail in clouds . Since the density outside the particles is much 
less than the density inside, we shall have to be more careful when considering emitted radiation 
from the surrounding gaseous medium. 

Under the assumption made, Equation 32 is acceptable as it stands upon a removal of the cav- 
ity walls. The temperature T will depend only upon the energy available through collisions, and 
the radiation field does not have to be in equilibrium with the surrounding medium. 

Averaging Techniques and the Equation of Transfer 

We are now in a position to evaluate the net change in intensity in the direction (/x, 4 ) due to 
the presence of the mass element dm . Adding up all the net gains and losses from Equations 16, 

21, 28, and 32 we obtain 



The terms on the right hand side of Equation 33 are negative or positive depending on whether they 
are respectively losses or gains. 
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After an examination of Equation 33, it becomes clear that a judicious choice of an effective 
phase function for single scattering p 0 (^, <j>; i±‘, <p‘) given by 


p 0 (m, 4>; . 4 >' ) 


/•CO /»CO 

Q s (l> 7,1 i 2 D i dl i = I Q s (i> itt ? D; p. O, <£; m' . <£' ) dr. 

J 0 ^0 


(34) 


will greatly simplify the equation. Upon multiplying both sides of Equation 34 by d os'/A-n, integrating 
over all solid angles, and remembering that £, p. <P‘, > 4 >' ) dw'/4n r = l, we obtain the result 


I 


p 0 (/*. </>; m' 


do/ 

v ) "4m _ 


(35) 


Thus, P 0 (m, <p; h', 4>‘) is normalized to unity just as p, (m, 4>; m', <£')• Upon inserting Equation 34 into 
Equation 33 and using the results of Equation 35 we obtain 


SI V (z,/x, 4>) 


N 0 dV 


do- I „ ( z - <£) 


/*00 

[Q s (i) + Q A (i) ] ^i 2 D i dr i 

J 0 


N 0 dv r f 

dcr QWvr i 2 D i dr i p 0 (m. <#>: 

•>0 •'o/ / 


da/ 


m' . 4 >‘ ) l„ (z. . 4 >‘ ) ~4^ 


N 0 dV 


*00 

+ B, (T) Q a < 1 > vr £ 2 D ; dr. . (36) 

*'0 

Upon using the relations Q E {i) = Q s (i) + Q A (i) , m = cos 0, and d« = sin 0 d0 d4>, Equation 36 becomes 


SI V (z, M. <£) 


N 0 dV 

d 5 ~ Jo Q E i> 77 r i 2D i dr i 


- I„ (Z, /J., <t>) 


477 


Jo Q / ) 7 Tr i 2 

Jo Q E (i> 7 rr i 2 


D. dr . 

i i 


D. dr . 

i i 



Po O’ & m' , <£')!„ (z, /x\ <£')<V d£' 


f 


Q s a > 7?r i 2 D i dr . 


J Q e ( i) 77r 1 2 ° 1 


B v (T) 


(37) 


Equation 37 may be made less bulky by defining an effective extinction cross-section Xe , an ef- 
fective scattering cross-section Xg , and incidentally an effective absorption cross-section x , 
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such that 


Xe 


X S 


X A 


f 

r 

Jo 

J »oo 

0 


Qjp) 77 r j 2 D i dr x 
Qgd) 77r ^ D. dr A 
Q. (l > 7rr 2 D . dr . 

A ill 




> 


The effective albedo for single scattering is then given by 


(38) 


cc 


0 


Xs 

x E 


(39) 


In the context of plane -parallel atmospheres it is clear that, for any layer of thickness dz, a 
volume element cylinder dv of height dz and base area da z is defined by 


dV - d cr z dz . (40) 

If now the linear dimensions of da z are made arbitrarily much larger than dz, but still small 
enough so that da z remains an element of surface area, then the sides of the cylinder may be neg- 
lected relative to do- z in determining the effective cross-section of dv when seen from a slant- 
path in the direction (a, 4). The cross-section do- of dv then becomes 

do - jj, d cr z . (41) 


Upon utilizing Equations 38 through 41, Equation 37 may be rewritten as 


Bl v (z, fj . , 4>) 

M -N oXE dz = 


4tt 



O, <£; 


n', <t>‘) I v (Z, n', &4>' - (l -» 0 ) B v (T) 


(42) 


It will be convenient to define a normal optical depth r v measured from the top of the cloud inward 
such that 


~ ~ NqXe^ ’ 


(43) 
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and a phase function for single scattering p(/x, <£; ^ , <£') normalized to & 0 such that 


pO, 4 >\ Po (m» <£; • 


( 44 ) 


Upon letting the ratio 8i v (z, /z, <f>)/dr v approach its limit as &r v - o, and upon replacing l v (z, /z, <£) 
everywhere with i(t, <£) and dropping the subscript v, Equation 42 becomes 


dl(r, /z, <£) 


n +i 

p(/x, 0; /z' , ) I(r, /z' , ) d/z f d<£' - (l-Sf 0 )B(T) 

l 


(45) 


It remains understood that Equation 45 is implicitly a function of the frequency of radiation v . 

Equation 45 is the equation of transfer for an arbitrary field of radiation. In general this ra- 
diation field is the diffuse radiation field which has originated from thermal emission from the 
cloud, surrounding atmosphere, and the ground. There will also be a contribution at relevant wave- 
lengths from the sun, which in fact is the sole source of radiation at visible and shorter wave- 
lengths. It will be convenient to distinguish between the reduced incident radiation 77F 0 e T ^ iM ° from 
the sun, which penetrates to the level r without suffering any scattering or absorption processes, 
and the diffuse radiation field which has arisen through one or more scattering and/or emission 
processes; here ttF 0 is the total flux in the frequency interval (z/, v +dv) of a beam of radiation 
crossing in unit time a surface of unit area normal to the beam, and /z 0 is the cosine of the zenith 
angle of the sun. The sun is assumed to be approximately a point source. 

It will turn out that the intensity of the diffuse radiation field is the quantity most easily han- 
dled in the equation of transfer, because the boundary conditions are much simpler to impose cor- 
rectly in this case than for the case in which the intensity of the total radiation field is the depend- 
ent variable. In order to separate the diffuse field from the directly transmitted radiation from 
the sun, we proceed in the manner given below. 

The Diffuse Radiation Field 

Consider a parallel beam of radiation from a point source at infinity of flux 77F 0 crossing a 
unit surface normal to the beam. The magnitude of the flux crossing a unit surface which is in the 
plane of the top of the cloud is 77/z 0 f 0 , where the cosine of the zenith angle of the point source is 
/z 0 , and this is (Reference 3) 


^ o F 0 



yUl(0, - /Z, 0) d<f) , 


(46) 
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where I (0, - /x, <p) is the downward intensity of radiation in the direction (-/x, <p). Since the only 
contribution is in the direction (- /x 0 , <fi 0 ) from the point source, the intensity i(0, - fx, 4 >) should be 
of the form 


1(0, -fi, <p) - f(M> s(m-^o) b(<£) Z(<f>-4> 0 ) . 


(47) 


where s(/x-^ 0 ) and S (<£-<£ 0 ) are Dirac delta functions, and f(/x) are g ( 4 ) are respectively functions 
yet to be determined of (/x) and (<£). Inserting Equation 47 into Equation 46 we obtain 


F 0 


f 


f*2jr 

M f (A 6 ) S (m-^o) d M J B(0) d< t> 


Mo f (Mo) g(<£ 0 ) 


From Equations 47 and 48 the downward intensities at the top of the cloud become 


(48) 


1 ( 0 , -Mt <£) “ S(^-/x 0 ) 8 (<^- 0 o ) . ( 49 ) 

The total intensity i(r, //, 0) associated with Equation 45 is the sum of the intensity i d (t, /x, <f>) 
arising from the diffuse radiation field and the intensity directly transmitted from the point source 
to the level r. By analogy with Equation 49 this latter intensity may be expressed in the form 

i t (t, m ,0> = ttf q 8(m +m 0 ) s (^"^o) h ( T ) . (50) 


where h(r) is a function of r alone yet to be determined. The total intensity then becomes 


!(r , /x, 0) - I D (r, fi, 0) + 77F 0 8(/x + m 0 ) 8(0- 0 q ) h(r) , 


(51) 


and the equation of transfer (Equation 45) after some reduction becomes 


dI Q (r, fi, 0) dhfr") 

> l dr + ^MFo S(/x +m 0 ) s [4> ~ < Po)~SF 2 = I D (t. M. 0) + wF 0 3(/x + Mo ) S(<£ -<£„) h(r) 


477 I P(M< 7>; M' . </>' ) I D (r, /x' , <£' ) d/x' dtf>' - ^ F 0 h(7 ) p(*x, d>; - /x 0 <£ 0 ) - (l -S 0 ) B(T) . (52) 

Jo J-i 

We suppose that dh(r)/dr and h(r) are in general non-zero. Then in the limit as /x - - /x 0 and -<^ 0 , 
Equation 52 becomes 


M—Mo [ 77 F 0 mS(/x + Mo ) s(d>-0 o ) dh i; ^ 

4^4> 0 


lim r 

o I^Fq 8 (/x +m 0 ) s (7> - 7 0 ) h(r) 

Tr'T'o 


(53) 
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yielding 


h( T ) = c 0 e" T/M ° , (54) 

where c 0 is the constant of integration. Upon replacing h(r) in Equation 50 with its equivalent 
from Equation 54 and letting r- 0, it is observed that (cf. Equation 49) 

c o s 1 ■ (55) 

In general then the equation of transfer Equation 52 for the diffuse radiation field becomes 


dl d (t,/x, 0) _ x 

M ” I D ( T > 40 ' 4^ p(m. </>; . <#>' ) I D (t, /*' . d<£' 

•'o J-l 

- T F o e ’ TAo p(^- ^ - (1 - 2 > 0 ) B(T> . (56) 

The contribution to the total intensity due to the directly transmitted flux may be obtained from 
(cf. Equations 46 and 50) 


F(t) 


F o e 


-y^o 


(57) 


It should be clear that Equation 52 could be modified to include n point sources by including n 
terms in the source function of the form 


- ~ f (i) 

4 r o 


-r/y.f 


p(i*- d; - M 0 (i) , d 0 (i) ) 


(58) 


where (i) refers to the i th point source. However, this shall be considered to be beyond the scope 
of this paper, even though in the far infrared terms of this form could be included to account for 
isolated emitting sources, such as discrete cumulus-type clouds. In any event Equation 56 will be 
the basic equation of transfer considered in this paper. Solutions to Equation 56 will be sought 
which describe the radiation field at the surfaces of plane- parallel clouds of arbitrary optical 
thickness which are bounded on both sides by a vacuum. These solutions will hereafter be refer- 
red to as solutions to the restricted problem. The solutions to the restricted problem will then be 
incorporated as functions in the formulation of certain integral equations which will describe the 
general problem; i.e., each vacuum will be replaced by a realistic atmospheric model which will 
include any ground or outside cloud effects, where relevant. This concludes the formulation of the 
basic equation of radiative transfer. 
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THE RESTRICTED PROBLEM 


Introduction 

The restricted problem deals with the evaluation of the angular distribution of radiation and 
net flux at the surfaces of a plane- parallel slab of particulate medium of arbitrary optical thick- 
ness which is bounded on both sides by a vacuum. There are two major tasks involved: 

1. the determination of the angular distribution of radiation diffusely reflected from and trans- 
mitted through the slab which is irradiated by an outside point source, and 

2. the determination of the angular distribution of radiation emitted by the slab from sources 
of radiation internal to the slab. 

It will be assumed that these two fields of radiation do not interfere and thus may be calculated 
independently of each other. 

The basic equation of radiative transfer which the total diffuse radiation field must obey is 
(Equation 56) 


dI(T, JLl, <f>) 


_ 1 _ 

I(t, fj. t 0) - 4^ 



p (^. 0; 


fJL 1 , <t>* ) I(r, jJL 1 , 4>‘ ) d/x' d f 


- 1 F 0 e T//i ° p(m, 4>; -/V <t>o) - ( x ~ z o) B(t) . (59) 

It is understood in Equation 59 that i(r, <£) refers to the diffuse radiation field; hence the sub- 
script D is dropped, as it will be throughout the rest of this paper. It should also be realized that 
all parameters in Equation 59 except the angular coordinates are dependent upon the frequency of 
radiation; therefore Equation 59 is valid only over a frequency interval in which all these param- 
eters are sensibly constant. And finally, the Planck function b(t) is replaced by b(t) ( v understood), 
which implies that the temperature T may be uniquely expressed as a function of r alone. 


The Equation of Transfer in the n fh Approximation 

The method of solution that will be followed throughout this section is the method of discrete 
ordinates in the n th approximation. Most of the work has already been well described by Chan- 
drasekhar (References 3 and 7), and we shall for the most part simply list the results and indicate 
the necessary equations which the reader would need in order to obtain numerical solutions. Some 
derivations are needed, and these are contained for the most part in Appendices A, B, and C. In 
particular, a derivation for the particular integral needed to complete the solution to the equation 
of transfer containing the Planck function is new, and this is given in Appendix B. 

In the method of discrete ordinates the continuous radiation field is approximated by 2n lin- 
early independent beams of radiation. Each of these beams is characterized by a unique direction 
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y i9 and there are a total of 2n of these directions, n each in the upward and downward directions. 
The discrete values of y k (i = ±1 , * • * , ±n) are picked such that the integral over y in the source 
function of the equation of transfer may M most accurately'' be evaluated by a quadrature formula 
containing 2n terms which is chosen for that purpose. The formula which will be considered in 
this paper is of the form 


( Mi ) % J f(M)d/x 


(60) 


where the s are the weights appropriate to a quadrature formula based on the division y L of the 
interval (-1, +1). This specific formula picked by Chandrasekhar is due to Gauss, where the 
values of y i are the roots of the even -degree Legendre polynomial P 2n (y), and the values of a i are 
determined from the expression 


a i - 

(V°) 

r P 2 n <M> 
1 M-Mi 

(i = ±1, •• 

1 ±n) 



J-l 1 




(61) 


The identities 


a i " a -i " "M-i 


(62) 


are valid for all values of n . 

This use of the Gaussian formula was criticized by Kourganoff (Reference 4) because it ap- 
proximates an integration over the whole interval (-1, +1), and the integral of interest (Equation 
59) is discontinuous at the boundaries of the cloud at y - 0. King and Florence (Reference 8) have 
recommended the Sykes double- Gauss method of fitting the Gaussian formula separately over the 
ranges (-1, 0) and (0, +1). The merits of this method appear to be quite sound and should be in- 
vestigated more thoroughly. 

Another approximation which will be incorporated consists of replacing the phase function 
p (y, 0; y' 9 0') determined from the Mie theory with a finite Legendre polynomial series expansion 
of the argument, cos 0, where 0 is the scattering angle; i.e., it is required that (Equation 4) 

N 

P (At, 0; m\0' ) = ?(COS0) % ^^^P^cos©) ’ (63) 

1=0 


where £ 0 is the albedo for single scattering and the values of 2o t (l = l, • • • , n) are constants in- 
dependent of 0 which provide a "best" fit to Equation 63 in the sense of least squares. A proof is 
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given in Appendix A showing that the series approximation in Equation 63 is normalized 
to 3 0 . 

From Equation 5 it is seen that 

P j ( cos ®) = Pj [^' + (l “M 2 ) 1/2 (1 -M /2 ) 172 COS (<£' . (64) 

Expanding Equation 64 in accordance with the addition theorem of spherical harmonics we obtain 


i 

Pj (cos 0) = P ; ( M )Pj (^') + 2 22 (TtSt Pf (AOP'OOcosmO*' 


(65) 


If it is now assumed that the intensity in Equation 59 may be written in the form 


N 

I (r, jj . , <f>) - I (0) (r, fj.) + y I< m > (r, fi) cos m ~ 4>) , 


( 66 ) 


it can be verified (Reference 7) by direct substitution that Equation 59 splits up into the (n + l) in- 
dependent equations 


dl<°> ( r, /x) 
^ dr 
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-T F o e " T/M ° 00 Pj (Mo) - (I-So)B(t) 

1=0 


(67) 


and 


M ^^ir’ ^ r IW(t, m) --§■ 7^ (TTT)! 2S »P"^)| i Pf (^)I W 

l = m 


N 

- T F o e‘ T/M ° (-l) I1ra ^TT§! SjPf 00 Pf (Mo) (m = 1, •••, N) . (68) 

l =m 


At this stage it will be expedient to define the form of the Planck function, b(t). It is not un- 
reasonable to suppose that in any realistic physical situation b(t) is sufficiently smooth and well- 
behaved so that it may be expressed exactly by an explicit infinite power series in r. A term-by- 
term differentiation of this power series would be expected to correspond exactly to the derivative 
of the original function. 
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In keeping with the spirit of approximating the phase function for single scattering by a series 
expansion containing a finite number of terms, N, we shall correspondingly approximate the Planck 
function by a power series also containing N terms. The function b(t) in Equation 67 is then re- 
placed by 


B(r ) 



(69) 


where the values of b r (r = o, • • * , n) are constants independent of r which provide the "best" fit 
to B (r) in the sense of least squares. 

In the method of discrete ordinates the integrals in Equations 67 and 68 are replaced by sums 
according to Gauss’s quadrature formula, and each of the (n + l) equations is replaced by an equiva- 
lent system of linear equations of order 2n . Solutions must be sought in approximations n such that 

4n - 1 > 2N . (70) 

The 2n (n + 1) linear inhomogeneous differential equations which replace the (n + l) linear in- 
homogeneous integro- differential equations given by Equations 67 and 68 are 


dlfl> (t, /x. ) 


dr 




z = o 






E 

L- J 
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(i - ±1, ***, +n; j - ±1, ***, ±n) (71) 


and 


^ ^ dr’ ^ "T^ (T^i 

l- m 


p r(^) I(m) ( T ’ M i) 


2 F o e 


N 

(TTm) I 3, P" ( Mi ) Pf (Mo) 

l=m 


(i - ±1, * * • , ±n; j - ±1 , ±n; m - 1, ■ * ■ , N). (72) 

The complete solution to the system of Equations 71 and 72 will involve 2n(N + 1) constants of in- 
tegration which can be determined from the 2n(N + l) boundary conditions 

I® ( 0, M-j ) = I <m ^ ( 0, M-i ) = o (73) 
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and 


I(0) ( r i’^i) = I(m) ( r i^i) = 0 (i = 1, n; m = 1, *■*, N) , (74) 

where r 1 is the normal optical thickness of the cloud. If the cloud is semi- infinite in extent the 
n(N + l) boundary conditions* 

e~ T i (TO (r, fi i ) -» e“ r I (m > (t, - 0 as t ^ oo ( i = 1 , * * * , n; m = 1 , • * * , N) (75) 

replace those of Equation 74 in order to secure the boundedness of the solution. 


The Complete Solution 

The complete solution to Equations 71 and 72 is required where the frequency of radiation is 
such that neither radiation from the sun nor thermal emission from the cloud can be neglected. 
Since both Equations 71 and 72 are linear, this solution may be obtained by first solving the homo- 
geneous parts of the equations, and then finding a particular integral, which, when added to the 
general solutions, satisfies the complete system of inhomogeneous equations. 


The Azimuth-Independent Solution. 

The solution to the homogeneous part of Equation 71 is (Reference 7) 


I <P) 
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( t ’M 


V" 1 M„ e"^ 


a= 1 


1 + ^i k a 
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i = 0 
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1 "Mi k a 


03 l ( k a) (Mi) 

1 = 0 


(i = ±1, * * • , ±n) , 


(76) 


*Often it is not clear that these conditions are sufficient, although they are always necessary. The actual criteria involved are that the 
formal solutions to the equation of transfer (Reference 3) 

f°° -rfiL. dr 

I(°. Mi) -J J(T, Mi) e — (i = 1, ••• , n) , 

** 0 'i 

where J(r, /x^ is the source function, do not become unbounded asT-®, For solutions in the n th approximation, J(r, /x^ contains 
terms of the form (Equations 71, 72, 92 and 98) M_ a e + aT and e aT which, when multiplied by (l/^.) e fX ' 1 and integrated over all 
t, do not become unbounded for all values of k a and k^ (Equations 79 and 104), as has been verified numerically by the author for 
several specific cases (& 0 5* 0). Thus the requirements 

M_ a = M-“ = 0 (a = 1, , n) 


needed to make all constants of integration determinate are not automatically insured by Equation 75* However, for all single scat- 
tering phase functions expressible as a finite series of Legendre polynomials, the form of the solutions are the same, and the require- 
ments M_ a = MJ^ = 0 (a - 1, * * * , n) shall be expected to hold for semi-infinite clouds in general. 
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where M ±a are the 2n constants of integration, and the quantities ( x ) are defined, apart from an 
arbitrary multiplying constant, by the relations 


(x) 



"A^A (*) d J,A ( x > 


(l = 0, •••, N) 


(77) 


where 


d i,a (x) 


j l 

2 1 + x 


(j = ±1, •••, ±n) . 


(78) 


The arbitrariness may be removed by defining (x) to be unity. Upon combining Equations 77 and 
78, and setting l = o, the characteristic equation for k results: 


l 


1 

2 



^(k)PA(Mi) 

+ Atj k 


(79) 


This equation is of order n in k 2 and admits, in general, 2n distinct nonvanishing roots which must 
occur in pairs as ± k a (a = 1 , * • • , n). It can be shown that the <f x ’ s obey the recursion relation 

21 + 1 - co-, 7 

= " k(i + l) it _ ITT £i-, , ( 8 °) 


where 


= 1 ; ^a ~ o (X < o) . 


(81) 


Equation 80 may be used to generate (l = 0, • • • , n) as a function of k, and Equation 79 subse- 
quently solved for the 2n roots ±k a (a = i, • • • , n ). In practice, for large n, it is probably more 
practical to solve Equations 79 and 80 interdependently on an electronic computer. 


In order to complete the solution, a particular integral must be found which, when added to 
Equation 76, satisfies Equation 71. In practice it will be easiest to find two particular integrals 
by setting first one, and then the other term of the inhomogeneous part of Equation 71 equal to zero. 
The sum of these two integrals thus found will then comprise the particular integral sought. 


N 

If the term in Equation 79 containing £ b r r r is suppressed, the resultant particular integral 
(Reference 7) is 


I. 


i (Mo) 


i_ y 0 e 
4 F ° 1 + /x./ Mo 
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(i - ±1, • • • , ±n) , 


(82) 
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where 


y 0 = h( Mo )h(- M o ) 


(83) 


ana 


n^o 


H(x) = 


H n. 

Q( 1 +k a x ) 


(84) 


Alternatively, if instead, the term containing e is suppressed, the resultant particular in- 


tegral may be expressed as 


MbCt)] = ^C r , s P s _ r ( Mi ) 


(85) 


where the values of C r s are constants (independent of r and ^.) which may be written as explicit 
functions of b r and 25 1 alone. The complete set of relations required to generate the various values 
of c r s (r ° o, • * * , n; s = r, ■ • • , n) are as follows: 


^N-l.N-l ^N-l 
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(r + 1 ) 
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(87) 

(88) 
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(89) 
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2(N - r ) - 3 


N-r-l 

2(N - r) - 1 


C r + 1 , N _i (r • 0, • • • , N - 2) 


(90) 


and 

C. . 


(r + 1) 


2(s - r) + 1 


(s-r + 1) (s - r ) 

2(s - r) + 3 C r+1 ,s+2 + 2(s - r) - 1 C r + l, S 


(r = 0, — ,N - 3; s-r+1, • • • , N - 2) 

(91) 


A derivation of Equations 86 through 91 and a systematic method of generating the various values 
of C r s are given in Appendix B. 


27 


The two integrals given by Equations 82 and 85, when added to the general solution in Equa- 
tion 76, comprise the complete solution to Equation 71, and this is 


I® (r, ^ ) 


1 +Mi k a 
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• MO 


(i - ±1 , * * • , ±n ) . (92) 

The 2n constants of integration M a (a = ±1, * * • , ±n) are to be determined in accordance with the 
2n boundary conditions 

I (Q) (0,M-i) = I (0 >(r I ,^ ; j = o (i ' 1, ' " * . n) . (93) 

In the case of conservative scattering (?„ r l) it is shown in Reference 7 that the solution to 
Equation 71 can be written in the form 
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(i - +1. * * • . ±n) , (94) 

where M a (a = ±1, • • • , ±[n - U), M 0 , and M n are the 2n constants of integration to be determined in 
accordance with Equation 93. The summations over a are restricted to (n - 1 ) terms, since Equa- 
tion 79 in this case admits only ( 2 n- 2 ) distinct and nonvanishing roots. 


If the cloud is semi- infinite in extent, the boundedness of the solution requires that the terms 
containing the constants M_ a in Equations 92 and 94 all be suppressed (cf. Equation 75), An applica- 
tion of the flux integral reveals further that M 0 in Equation 94 becomes 4miSy in this case, as is 
demonstrated in Appendix C. The solutions contained in Equations 92 and 94 in the semi-infinite 
case then become respectively 
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and 


M a e a 

1 + Mi k a 


l( +k a)Pj(^i) + l+ Mi /Mo ^ 




3 ^o [l 1 " 3 S i) T + Mi + M n ] J' ( S o 


= 1; i = ±1, •••, ±n) . (96) 


The constants M a (a = 1, • • • , n) in Equation 95, and M a (a = 1 , • • • , n - l) and M n in Equation 96, 
are to be determined in accordance with the n boundary conditions 

l(°>(o, = 0 ( i - 1 , * * * , n) . (97) 

The Azimuth- Dependent Solutions. 

In a fashion entirely analogous to that of obtaining the solution of Equation 71 it can be shown 
(Reference 7) that the complete solutions to the azimuth-dependent equations (Equations 72) can be 
written in the form 

r * - k ™T 
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•— £ =m 

where the 2nN constants M ± * are to be found from the conditions 

l( m ) (o, fi _ i ) = (tj , fi i ) = 0 (i - 1, m - 1, ***,N) , (99) 

and r l is again the normal optical thickness of the cloud. If the cloud is semi- infinite in extent, 
the boundedness of the solutions requires that the terms in Equation 98 containing M™ a all be sup- 
pressed. In this case the solutions become 
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where the nN constants M a m are found from the conditions 


I (m) (0, M-i) = 0 (i = 1, •••, n; m = 1, •••, N) .. (101) 

The quantities .f j" (x) are defined, apart from an arbitrary multiplying constant, by the relations 


^f(x) 
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The arbitrariness is removed upon defining £ m m (x) to be unity. The resulting characteristic equa- 
tion, upon setting l = m, becomes 
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£>A (X~+ m j ! ^ M P ra m (Mj) P A m (Mj) 
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1 + Mi k” 


(j = ±1, •••, ±n) . (104) 


Equation 104 is of order n in (k m ) 2 and admits, in general, 2n distinct nonvanishing roots which 
must occur in pairs as 


±k a m (a “ 1 , • • • t n; m = 1 , • • • , N) . 

It may also be verified that the <f A m ’s obey the recursion relation 

2I + 1 "“> , Z + m . 

£ m = — £ m _ ~ £ m 

s;+l k m (Z-m + l) ?i Z-m+lM-i 

where 


(105) 


C = 1 ; ^ = 0 (A < m) . (106) 

As before, Equation 105 may be used to generate the functions of £ A m (k m ) required for the evalua- 
tion of the roots of Equation 104. Once these roots have been obtained, the values of <f A m (k m ) follow 
from Equation 105. 
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The values of y 0 m in Equation 98 are made determinate from the relations 


r 0 m = H<*> (Mo) H (m > (-Mo) (m = 1, ' ’ • . N) , 


where 


H< ra > ( x ) 


n ( x+ ^i) 

i=l _ 

11(1 + k"x) 

a - 1 


( 107 ) 


(108) 


The Angular Distribution of Outgoing Radiation 

It is required to evaluate Equation 66 for the angular distribution of outgoing radiation at the 
surfaces (r = 0 and r r tJ of the cloud, and this may be accomplished through nothing more than 
simple substitution. The form of the solution will depend upon the number and type of sources of 
radiation involved. It will be convenient in the present context to treat each source separately and 
add the separate resultant solutions to obtain the final solution. This will entail no loss in gen- 
erality and will circumvent any ambiguities which might arise later in the treatment of the general 
problem (in the next section), where the bounding vacuums in the restricted problem are replaced 
by more realistic physical conditions. The only assumption required is that the separate radia- 
tion fields do not interfere with each other. 

In order to reduce the space required to write the various solutions, the following artifices 
will be employed: 

1. The azimuth- independent and dependent terms will be combined by extending the summa- 
tions over m to include (m - 0, • • • , n); e.g., in this abbreviated form Equation 66 becomes, 
in the scheme of the n th approximation, 
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■ r I (m) (t, cos m(0 o - <£) ( i - ±1, • • * , ±n) . (109) 

m = 0 

2. The symbols 8^ ^ and Sj" will be respectively defined by 


8 


/3.\ 




(\ = /}> 
/ & 


} 


and 


( 110 ) 


~ Ll ~ m ) ! 

w l ~ a i (1 + m) ! ' 


( 111 ) 


31 



3. The definition 


k ® 


(a = 1 , 


n ) 


(112) 


will be employed, and the symbol ^2 j will be understood to mean the summation over a, 
where the range will be specified in accordance with the relevant solution. 

With these definitions in mind we now turn to a description of the angular distribution of radiation 
at the boundaries of a plane- parallel cloud within the context of the restricted problem. 


Diffuse Reflection and Transmission . 


The solutions relevant to radiation from an outside point source which is diffusely reflected 
from a cloud of normal optical thickness are (cf. Equations 92, 94, 98, and 109 through 112), 
for a 0 / 1 :• 
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If the cloud is semi- infinite in extent, the range of summation over a in Equations 113 and 114 must 
be restricted to (a = 1 , - • • , n), and the constant M 0 in Equation 114 identified with unity (cf. Equa- 
tions 95 and 96). In all cases the usual relevant boundary conditions (Equations 99 and 101) have 
been imposed in determining the various constants of integration. 

The corresponding solutions for the radiation diffusely transmitted through a cloud of normal 
optical thickness r l through analogous considerations become respectively, for 2> 0 / l; 
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and, for S 0 - 
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Thermal Emission and The Law of Darkening . 

In view of the fact that no outside sources contribute to the diffuse radiation field in our pres- 
ent context, the solution for the angular distribution of outgoing radiation from a thermally emitting 
cloud reduces to an axially symmetric one, and the azimuth- dependent terms in Equation 109 are 
all suppressed. From Equation 92 it is seen that the solutions at the top and bottom of a cloud of 
normal optical thickness t 1 are, respectively, for r - 0: 
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and, for r - r, : 
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If the cloud is semi-infinite in extent, it can be seen from Equation 95 that the range of integra- 
tion over a in Equation 117 must be restricted to (a = 1, * * * , n). 

It is clear from Equation 67 that in the special case of conservative scattering (a> 0 = l) in a 
semi-infinite cloud the law of darkening is independent of the temperature. Since no sources or 
sinks are available to the radiation field the net flux is constant (independent of r). In terms of 
this constant net flux rrF the solution to Equation 71, upon suppressing the inhomogeneous terms, 
can be shown to be (Reference 3) 
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from which it follows that the law of darkening is given by 
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The values of ^ and k a may be determined from Equations 79 and 80 and the constants 
M a (a = 1, * • • , n - 1) and M n from the n boundary conditions 

l(0, - M .) = o (i = 1, • • • , n) . (121) 

It should be noted that any time- independent solution for the angular distribution of conservatively 
scattered radiation in a cloud of finite optical thickness is meaningless. 

Molecular Effects 

In general it may be supposed that the atmosphere within the cloud contributes to the absorp- 
tion and thermal emission of radiation, but not significantly to the scattering. Any mass element 
dm in the cloud then absorbs radiation through both particle absorption and absorption by the free 
molecules contained in dm. 

Let E be the radiant energy incident on dm from all directions in unit time which is extinguished 
(absorbed plus scattered) by dm. If g is that fraction of E which is extinguished by the particles 
contained in dm, then (1 - g) must be that fraction of E which is absorbed by the free molecules in 
dm. If E p is the amount of energy extinguished by the particles in unit time, then 

Ep s = a 0 E p = S 0 gE (122) 

is the amount of energy scattered by these particles in unit time, and 

V = ( 1 _ 2 ; o) e p r (l-" 0 )gE (123) 

is the amount of energy absorbed by the particles in unit time. It is further clear that 


E m A = (1-8)E 


(124) 


must be the energy which is absorbed by the free molecules in dm in unit time. Upon adding Equa- 
tions 123 and 124 it is seen that the total amount of energy absorbed by dm in unit time must be 


V + e m a = (i-S3 0 g) E - 


(125) 


The ratio of the energy scattered in unit time to the energy extinguished in unit time, i.e., 
the effective albedo for single scattering of dm, is given by (Equation 122) 
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If further, local thermodynamic equilibrium is assumed to prevail in dm, the rate (5) of emitted 
energy in accordance with Kirchhoff's law is seen from Equation 125 to be 


3 = (1-S3 oB )b(t) • (127) 

And finally, since the effective extinction cross-section of dm has been increased by a factor of 
(l/g), the differential optical thickness dr becomes (Equation 43) 


dr 


~ I N 0 *E dz 


(128) 


If g is independent of r, it should be clear from Equations 42 and 126 through 128 that, upon includ- 
ing the effects of free molecules upon the radiation field, a system of solutions results which is 
equivalent in every respect to the system of solutions discussed previously, merely by replacing 
y e and a) l ( l = 0, . . . , n) everywhere with (l/g) y E and 2 j z g ( i = o, - • - , n) respectively. 


THE GENERAL PROBLEM 
Introduction 

In order to assess properly the radiation leaving a planetary atmosphere, the interplay of 
radiation with the gaseous atmosphere and aerosols, as well as with clouds, must be taken into 
account except perhaps for limited regions of the electromagnetic spectrum. To further compli- 
cate the picture ground effects must be considered where clouds are optically thin. 

Several points of view are possible in any treatment of a physical problem. Ours shall be 
one directed toward (though not restricted to) practical numerical solutions connected with certain 
physically realistic problems associated with the transfer of radiation through the Earth's cloudy 
atmosphere. Since our interest is primarily utilitarian in nature, it will be necessary to be some- 
what restrictive in the formulation and solution of the problems considered. In particular the nine 
i approximations enumerated in the Introduction will be maintained. 

Two basic problems are considered in this section. The first problem considered is the de- 
termination of the intensity distribution of short wavelength radiation diffusely reflected from an 
j optically thin Rayleigh scattering atmosphere overlying an optically thick cloud. The second 

problem is the determination of the intensity distribution of long wavelength radiation at the top 
j of an atmosphere containing an optically thin cloud. In the first problem an outside point source 

(e.g., the sun) is assumed to illuminate the atmosphere- cloud composite from above. In the 
second problem the radiation is assumed to arise from the atmosphere- cloud- ground composite 
i-V by thermal emission associated with a characteristic temperature T at an optical depth r. in the 
latter problem the ground is assumed to emit radiation in accordance with its characteristic 
blackbody temperature T. In neither case is a contribution due to aerosols considered, although 
T: an obvious extension of the problem could be carried out to include them. In view of the utilitarian 
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purpose of this paper it is felt that further extensions of this sort (including optically thin clouds 
at short wavelengths) would lead us too far from the main theme, since the numerical solutions 
would rapidly become overwhelmingly complex. 

Diffuse Reflection and Transmission 

Consider a plane- parallel slab of a gaseous or particulate medium of normal optical thickness 
r 1 bounded on both sides by a vacuum. We shall want to define a scattering function 


and a transmission function 


s ( T r M- <t>\ M 0 - <t> 0 ) 


l( T i > M- <t>\ M 0 - <£ 0 ) 


such that the diffusely reflected intensity of radiation at the top of the slab is given by 

lo , 

1(0, M- 0) - 4/i S ( T 1- M, 4>' Mo- <P 0 ) - (129) 

and the intensity of radiation diffusely transmitted through the slab at the bottom of the slab is 
given by 


0 

I ( T l* -M, <*>) = 4^ T ( t i ; <P< Mo- <t>o) ■ (130) 

These intensities are solutions to the so-called restricted problem. The radiation in Equations 
129 and 130 is assumed to be originally in the direction (~av ^o) and thence redirected respec- 
tively into the directions (m, <P) and (-m,<£ ), where the cosine of the zenith angle m is in the range 
(0 <m< 1 ) and the azimuthal angle 0 is in the range ( 0 < <P< 2n). The quantity F 0 is related to the 
net flux ttF 0 normal to the incident beam of radiation. 

It should be noted that the reflected and transmitted intensities refer only to the radiation 
which has suffered one or more scattering processes. For example, if we denote the direction of 
incident radiation by ^o)> I ( T i’ $) does not include the directly transmitted flux ttF 0 e Tl /^° . 

The particular forms of the scattering and transmission functions in Equations 129 and 130 
were chosen from the following considerations: 

If the intensity of an outside point source is given by Equation 49, it follows from the flux in- 
tegral (Equation 8) that the net flux diffusely reflected by the slab is 





\ d&j 

M, <P; Mo- 4> 0 ) 4tt 


(131) 
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and the net flux diffusely transmitted through the slab is 

T(r t ; /*„, 4> 0 ) % , (132) 

where, from the boundary conditions, 

s ( r i ; 0; M O ’0 O ) = T ( r i’ M 0 » ^o) = 0 (0</x<l) , (133) 

* and the integrations are to be performed over all solid angles. Thus S (r x ; /j, 0 , 0 O ) and 

T (rj; <p; jj. 0 , 0 O ) are normalized so that their integrals are respectively a measure of the fluxes 
diffusely reflected from and transmitted through the slab relative to the flux of incoming radia- 
tion crossing a unit area normal to the incoming beam. 

The scattering and transmission functions as defined obey Helmholtz's principle of recipro- 
city; i.e., 

S(t x ; 0; m'. 0') = s(r x ; 0) (134) 

and 

T(r x ; (j., <f>; /jl' , <f>*) = T (t x ; /z\ 0' ; jj,, 4>) , (135) 

where the more general primed quantities replace the explicit values and 0 O . With these defi- 
nitions in mind we turn now to the formulation and solution of the integral relations which express 
the intensity of diffusely reflected radiation at the top of the atmosphere in terms of the restricted 
problems. 

Formulation of the General Problem . 

i 

The solutions to the restricted problems suggested by Equations 129 and 130 will be presumed 
to be known. For example, Chandrasekhar (Reference 3) has given the restricted solutions to the 
! problem of diffuse reflection and transmission for a plane -parallel slab of gaseous atmosphere, 

* and the comparable solutions for a plane-parallel cloud layer have been indicated in the previous 
section of this paper. 

> 

The problem now is to unite the solutions of the restricted problems with the scattering and 
transmission functions which arise as formal solutions to the "composite" problem in such a way 
as to allow a tractable solution for these functions. It will turn out that two interdependent linear 
■ ; integral equations plus one linear algebraic equation suffice to formulate the problem. Figure 3 

1 0 illustrates the role of each intensity component which is considered in the following discussion. 
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Let 



Ia( 


°-v: 


,! W 


1 


I A (0, M. <t>) = V S A (r,; fi, 4>\ M 0 . <*> 0 ) 


(136) 


and 


Figure 3 — The intensity components arising through the F 0 

interaction between an incident beam of parallel radio- I A (r t , -/x, 0) = T A (tj ; /jl, 0; jj, 0 , 0 O ) (137) 

tion of flux 77f 0 and a model of the atmosphere -cloud 
composite. The tops of the cloud and atmosphere are 
situated respectively at r = r 1 and r = 0. 

be the restricted solutions for the intensities at 
the top and bottom of a slab of gas respectively 
as derived by Chandrasekhar (Reference 3). 

Further let 


* 0 . . 
I c (0, /x, 0) - 4/2 S c (/X, 0; /x 0> 0 O ) 


(138) 


be the restricted solution for the intensity at the top of a semi- infinitely thick cloud due to diffuse 
reflection from the cloud. 

In what follows the subscript "A" shall always refer to the atmosphere above the cloud and 
the subscript "C" for the cloud itself. All quantities with an asterisk will refer to unknown in- 
tensities arising from the interplay of radiation between the cloud and atmosphere. All quantities 
without an asterisk will refer to known intensities arising as solutions to the "restricted" problems. 
All quantities with (~) above them indicate a reversal of direction, or a so-called change of parity; 
for example, if T denotes the transmission function through a slab in one direction, T denotes the 
transmission function in the opposite direction. 

Let I* (tj , }jl, 0) be the intensity in the direction (^x, 0) at the top of the cloud arising from the 
interaction of radiation with both the atmosphere and cloud, and l A * 0)be the corresponding 

intensity at the top of the cloud in the direction (->, 0 ). The normal optical thickness r 1 of the at- 
mosphere is the optical depth at which these intensities are located. 

The following equations are somewhat more general than those of Chandrasekhar, where he 
considers the more restrictive case of an isotropically scattering ground in his discussion of the 
"planetary problem" (Reference 3, page 269). The derivations of the formulas parallel his. 
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I* ( r i» At. 0) will be composed of two parts. One part is due to the radiation directly trans- 
mitted through the atmosphere before being diffusely reflected by the cloud. Symbolically we may 
indicate this component as 


I, = e'V M « I c (0. M. <f>) . (139) 

where the multiplying factor e Tl/ ^° is due to the attenuation the directly transmitted radiation un- 
dergoes in traversing the atmosphere of normal optical thickness r x . 

The second component is due to the intensity of radiation from any direction in the atmosphere 
incident on the cloud top which is diffusely reflected by the cloud into the direction fa, 0 ). 
Symbolically 


I 


2 


4tt JT f f S c O, 0; ft', 0 f ) X A * , -V , 0 # ) <V d0 # . 


(140) 


The integration is performed over all solid angles in the upward hemisphere in order to ac- 
count for the radiation from the sky incident on the cloud from all directions. Note that I*^, -//, 0') 
is the ’’unknown” intensity of sky radiation at the cloud top in the direction (-/ 0'). 

The intensity l A * (t 1? 0 ) is similarly composed of two parts. One component is due to the 

intensity of radiation diffusely transmitted through the atmosphere in the direction (->, 0 ) which 
originates with the main source (e.g., the sun). Symbolically we have 

F o 

1 3 = 4 >) = t a Cv- 4>\ m 0 . 4>o) ’ (141) 


the transmitted intensity derived from the restricted problem. 

1 The second component is due to the intensity of radiation (r 1( p , 0') in the direction (/x',0') 

arising from the cloud top, which is diffusely reflected by the atmosphere into the direction (-7X, 0). 
) Symbolically we have 



-/X, 0; -{1 * , <£') I c * ^T x , fx * , 0') d \x d0 f 


where the integration is carried out over all solid angles in the upward hemisphere. 


(142) 
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Upon collecting Equations 139 through 142, we obtain two independent integral equations relating 
the intensities l c * (tj, /x, <£) and I A * (r x , -p, <£) at the top of the cloud, namely 


IcK.m.*) = 


1 0 


lA, 


j_ rr 

Jo •'O 


s c(^-^ M o . 0 o ) + 4^1 s c ( M ’ <t>i ( 7 


0') d/x' d0' (143) 


and 


( T i ’ "M. 4>) 


4 M t a( t i'- 0 : /v 0 O ) 


n l 

S A (r i; -M. 0: V. 0') Ic (-1 . m'. 0') V «#' .(144) 


Upon inserting Equation 144 into 143 we obtain 


/•277 /%1 


I c ( 7 i. M. 0) 


^0 -T, /M n / \ 1 

7 s c (m, 0; M 0 . 0 O ) + 4ttm 


S c (m, 0; /z\ 0')f^T r T a(t 1: ^'.0'; Mo- 0o) 


A 7 "! J S A (r i; -m'. 0 '; -m". 0") I c * ( T l- /*"• 0") dM"d0*lcV' d0' . (145) 


Equation 145 is an integral equation for I c *(t 2 , 0) expressed in terms of the known scattering 

and transmission functions of the restricted solutions. 

Finally, we must obtain the intensity of radiation in the direction (^, <£) at the top of the at- 
mosphere due to the presence of both the atmosphere and cloud. This intensity l*(0, ^,<p) is com- 
posed of three parts. The first part is the radiation from the source diffusely reflected by the 
atmosphere, I A (0, <fi), as given by the restricted solution. The second part is the intensity 

arising from the cloud top in the direction (^, <£ ) directly transmitted through the atmosphere and 
attenuated by the amount e Tl/A \ Symbolically this intensity is 


= I r 


( r i- 4>) 


l/M 


(146) 


The third component is the intensity of radiation arising from the cloud top in the direction , 4>') 
and diffusely transmitted through the atmosphere into the direction (m, <£). Symbolically this 
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intensity is 


1 6 " ~^ji\ T a (r,; , 4>‘) 1* (t , 4>‘) d^' Af}' . (147) 

‘'0 •'0 

Once Equation 145 has been solved for 1^ /x, 0), the outgoing intensity I* (0, /x, 0) at the top of 

the atmosphere may be found with the aid of Equations 146 and 147 and the related discussion. Thus 

I* (0, /x, 0) = I A (0, /X, </>) + I c * (t x , /X, 0) e 1//M 

J* 2 tt /• 1 

T a ( T ii “M. <£: ~M*. 0') I c * (t,, m' . </>') <V d0' . (148) 

0 *'0 

Equations 145 and 148 represent the formal solution to the outgoing intensity of diffusely reflected 
visible radiation in the direction (/x, 0 ) at the top of an optically thin atmosphere overlying an 
optically thick cloud. 


Solution to the Problem . 

The formulation of the problem is contained in Equations 145 and 148. In keeping with the 
definition of scattering functions we may define a scattering function such that (Equation 129) 

F q 

I C *( T 1 ,M. 0) = s c * (Tj! m, 4>\ /v 4>o) ■ (149) 

Upon inserting Equation 149 into 145 we obtain 


S c ( T r /a 0; 


n l 

S c (m.0; m'. <t>') ~t r T a ( t l ; M*. <*>': Mo-^o) d ^' d<£' 


/*2rr/»l 

7717 I S c (^,<#>; /,(#>') 4 S A (r i; V,0'; Mq.^o) df/’d^'d/i' d</i' . ( 15 °) 

*^0 0 ) 0 ) •'0 


We shall require the solution for S* 0; /u 0 , 0 O ) from Equation 150. 

The restricted solutions in Equation 150 have been solved in general in the so-called n th ap- 
proximation in the form (Equation 109) 

N 

r ( T ' <■.•*) = r I (m) (r, jLt cos m (0 O - 0) ( i - ±1, • • * , ±n) , (151) 

m — 0 
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where fx ±i are the intervals appropriate to an n-point quadrature formula picked for the purpose 
of evaluating the integrals which appear in the source function of the equation of transfer. 

In keeping with Equation 151 we re-write the relevant scattering and transmission functions 
in the form 


and 


S c ( M , <*>; m' . <*>' ) = 2^ s c (k) o- m') cosk(0' -<P) , (152) 

k = 0 

N 

S A ( t j; -/a', 4>'\ ~m". 4>") - £V><V.-m')cosI ( *W') . (153) 

1 = 0 

N 

s c * (Tji /A 0 , <£ 0 ) = 2^ s c* (m) (m. Mo) c°sm(<£ 0 - 0) , (154) 

m= 0 


N 

T A ( t i; m*. 0'; m 0 - ^o) = £>"> (m'.Mo) cpsn K-0 . (155) 

n = 0 


where the dependence on r x is not indicated merely for the sake of convenience. 

Upon replacing the second term on the right-hand side of Equation 150 with we have suc- 
cessively with the aid of Equations 152 and 155 the relations 


■ 277 /»1 

s c (m> 4>\ t a ( t j; m'.<£'; m 0 > 0 o ) d M* d< £' 

•'O ■'0 


/•2 77 

• 1 

N 


N 



r s c k) (^’ ^ c ° sk (^ # ~4 > ) 

1 

i 

f 1 

71 T A <n) (m'. M 0 ) cosn(0 Q -0') 

Jo « 

0 

__k = 0 


_n = 0 


dpt' d0' 



r 27T 

C 1 

N N 

JL_ 

477 

0 • 

JL 

0 1 

22 22 Sc<k) ( ' m ’ m *' > Ta(ii) cos k ^' cos n ^° 

_k=0 n = 0 _ 


d^ f d <£> ' 


(156) 


With the aid of the relation 


r 277- 

COS k(<£' - cp) cos n(<£ 0 - <£') d <j> 9 

Jo 


fo (k / n) 

<77- cos n A £ 0 - (k - n f 0 ) > 

(j2n- (k = n = 0)J 


(157) 
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Equation 156 becomes successively 


/•i 


J7 

477 


1 


zz s c < k > (M, M')T A ( n >( M ', “"'I cos k(0' ~ 0) cos n (0 Q - 0' ) d0' 


d/x # 



(l +S o.n) S c <n) (M. m')T a <")(m',Mo) 


d/7 


cosn(0 Q -0) , 


(158) 


where 


^ 0, A. 


{» 


(\ = 0) 
(\ / 0) 


(159) 


In like manner, for the third term on the right-hand side of Equation 150, denoted by i 2 , we 
obtain with the aid of Equations 152 through 154 successively the relations 


r2-n fl f2rr fl 


1677 2 


-4 

1677 2 I 

•Jq 


S c (m. <t>\ /A 0') 77 S A (r x ; -/a\ 0'; V' . 0")t 7 S c ( r i ; M". <£"; <£ 0 ) d <£" V dc £' 


N ”1 N ^ 

'Y^ t S c (k) (At. A 1 ') cos k(0 # ”0) 77 YY S A (i) (”/*'. *V') cos *(<£" ”<£') 


L i = o 


b •'o •'o •'0 


is 

7 S c (m> (^** M 0 ) cos m (cfi g - <p"J 


d/i" d^6" d/i’ dtp' 


r 2 n r‘ r 2n 1 

id Njf*[££ 

Jo Jo c° 0 Lk = 0 1 = 0 


S c (k > (fJ., fJL* ) S A < ( -fl* , ) cos k(0' - 0) COS £(0" - 0 1 ) 


d/7 d0 # 


r N 


7 S c* <m> (/*"• ^ 0 ) cos m ($0 


dfi" d 0" . (160) 
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By analogy with the derivation of Equation 158, Equation 160 becomes successively 



x cos Z(0" - 0) cos m (0 Q - 0") 


C V d0" . 


d 

t 

V- 



Upon replacing k, i, and m by n and utilizing Equations 152, 154, 158, and 161, Equation 150 
becomes 


N 



n = 0 


cos n(0 o -0) 


N 



cos n (0 O -0) 


n = 0 



>(^0 ~$) 


S c (n > (/*. m') T a <"> (m'. Mo) 


d jjL t 



-*) 


s c *C"> 


(m" 



s c (n > (M, /x')S( n )(- M ', ~/i") 



(162) 


The exponential in the first term on the right-hand side of Equation 162 serves to remind us that 
each scattering and transmission function which is subscripted with an TT A" or superscripted with 
an asterisk is an implicit function of the optical thickness r 1 of the overlying atmosphere. Let us 
define 


f 


Sj") (fjL, m')T a <"> (m' 


Mo) 




- h/ n >( M , M 0 ) 


(163) 
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a ’’known” function of (^, /x 0 ) for a given r l , and 


I' 


S c (n > (m. M*) S A <") (V. -M*) ^ 

^ fl 


= hC n > (At , M ") 


a ’’known” function of (a^m" ) for a given r 1 . 

With the aid of Equations 163 and 164, Equation 162 may be re-written as 


N 

^ ' s c* (n) (m. M 0 ) cosn (^> 0 - 0 ) 


r i/ M o 


s c (n) (m. Mo) cos n (<£ 0 -<p) 


N 

(l + 5 0 .n) h i (n) ( M ’ ^o) cos n (<^> 0 - 


*) 


N 


> (M. m") S*<») (m". Mo) 


cV 




Equation 165 must be valid for all 4>. We therefore have a set of (N + l) equations of the form 

S c * ( "> (r.i M, Mo) = e’V-o s< (n> ^ ^ + | (l +8q J h/") (r,; M, M 0 ) 

+ H( 1 +S o , n ) 2 f M(M '')S*( n %i; m".M 0 ) ^ (n = 0, ••• 

•'O 


N) 


where the dependence of the various terms upon the parameter r x has been re- indicated. 
Equations 166 may be written in the more concise notation 


S c * (n) (T,; Mo) = f (n) (t,; m. M 0 ) + f g (n) ( T i; M. m') S c <n) ( r 1 ; ^ ’ M 0 ) d M' ( n ~ °> 


. N) 


where 


f< n > (Tj; m. M 0 ) - e Tl / M ° S c < n ) (m. M 0 ) + T l 1 + 8 0,n) h l (n) ( T l : M o) 

g< n > (TjJ M, M') = T^I 7 (l +S 0,n) 2 h 2 <n> ( T l= M. M') • 


and Tj and At 0 are given parameters, constant for a given problem. 


(164) 


(165) 

(166) 

(167) 

(168) 
(169) 
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Equations 167 are formally examples of Fredholm's linear integral equation of the second kind 
and may be solved by standard methods. One method which lends itself well to the particular 
problem at hand is presented in Appendix D. 

The final quantity which we wish to solve for is the outgoing intensity in the direction (m, 0) at 
the top of the atmosphere given by (Equation 148) 

I* (0 , /x» 0) = I A (0, /x, <p) + e I c * ( T i . 0) 


n i 

T a ( t i; ->,<*>; V. 0 ') I * (tj, /x', 0 ') d/x' d 0 ' 

J 


or, in terms of the relevant scattering and transmission functions, 

S* ( T i » <p; M 0 1 <t> o) r S A (Tj; m. o) + e TlA S c * (r,; M, 0; M 0 . 4> 0 ) 


n l 


T'a ( t i ’ “A 6 * 0! S c ( T i ; M 0 . ^o) d0' , (170) 


where 


* o * 

I (0, /x, 0) = 4^7 S (t,; m. Mo.'/’o) ■ (171) 

Upon denoting the third term on the right-hand side of Equation 170 by £l 3 , we find with the aid of 
Equation 154 and the relation 


iv 

T a (ti ; -/X, 0; -/J.' , <p') r ^ T A (n > (“/X, "/X* ) cos n(0' - 0) 


that ^ 3 may be successively written as 


477 


t A ( T i ; '»<£') 77 s c * ( t i^ a^'» 0 '; av <£ 0 ) d M' dc £' 


(172) 


477 

- 2 tt 

r 1 

N 

V -1 T A (n) ( - fjL ' ) cos n (0' -0) 

1 

‘m 7 ’ 

N 

S*<-> (m' > Mo ) cos m(0 o " <£') 

J 

0 * 

0 

n = 0 _ 


m= 0 
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(*2 7T 


JL_ 

4 tt 


T A (n) ( "m') s c* (m) (m'* Mo) cos n (0' ”0) cos m (0o “0') 

_ n= 0 m=0 

By analogy with the comparison of Equations 156 and 158, we obtain successively 


djj,' d0' 


(173) 


N 

2] ( 1+8 0.n) T A (n) (-M. V)S*("> (m',M 0 ) cosn(^ o -0) 


d/x f 

M # 


(174) 


or 


where 


0(n) 


(m- Mo) 


N 

3 3 = 0< n) (m. Mq) cos n (0 O - 0 ) 

n = 0 


T( 1 + s o.„)J o t a <">(-m, _ m’> s c * (n) (m\ Mo) 


(175) 


d^i* 

7^ 


(176) 


From the manner in which it enters Equation 170, we are led to assume that the scattering function 
s * ( r i5 M, 0; m 0 ,0 O ) ma Y written in the form 


IN 

s * ( r i ; 0; 0 O ) = 22 S* (n) (M* /^o) cos n (^o • 


(177) 


n - 0 


Upon gathering terms from Equations 170 and 175 we find with the aid of Equations 154 and 177 
and the relation 


S A (t i; m , < P - Mo , 4>o) - S A (n) (M. M 0 ) cosn(<^ 0 - < f >) 

that Equation 170 may be written in the form 

N N 

22 s * (n) M 0 ) cos n ( 0 Q - 0 ) = 2 ~ 1 | S A (n) (m. Mo) cos n( 0 Q - 0 ) 


(178) 


+ e 


™ N 

22 S c* (n) (m* M 0 ) cos n (0 o “0) + 2~! £> (n) (m* Mo) cos n (0Q -< £) * (179) 
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Since Equation 179 must be valid for all 4 > , it separates into (N + l) equations of the form 


s*<"> (T,; m. Mo) = S A <"> (rj ; /*, m 0 ) + e‘ 


lA 


S*<"> 


M 0 ) + ® (n) ( T j; M. Mo) ( n " 0, 


N )> (180) 


where the dependence of the various terms on r 1 has been re-indicated. A method in keeping with 
the context of the problem for evaluating the integral 0 (n) (rj /x, /x 0 ) in Equation 176 is given in 
Appendix D. 

One remark should be made at this point. If the solutions to the restricted problems have 
been determined in the n th approximation, the solutions to the composite problem as indicated by 
Equations 167 and 180, for example, should also be solved in the n th approximation (cf. Appendix 
D). This completes our discussion of the problem of diffuse reflection of visible radiation in an 
optically thick cloud plus overlying atmosphere composite. 


The Law of Darkening 

The second problem which we shall investigate is the problem of limb darkening in the in- 
frared region of the spectrum at wavelengths of a few microns where an optically thin cloud is 
present. There is nothing in principle that makes the solution for the outgoing intensity more dif- 
ficult for a finitely thick cloud as far as the restricted problem is concerned. One simply re- 
places the n boundary conditions (Equation 75) 


l(r, /x i , 0je T -»O as r 00 


(i - 1, * • • , n) 


with the n boundary conditions 


I ( T i- Mi. 0) - 0 (i - 1, • • • , n) , 

where r x is the normal optical thickness of the cloud. It is now necessary to solve 2n linear al- 
gebraic equations for 2n unknowns in this case. 

The problem becomes considerably more bulky when ground and atmospheric effects are in- 
corporated into the results. However, if one is willing to admit the approximation that the ground 
emits as a blackbody at its characteristic temperature and that aerosols and dust in the atmos- 
phere contribute a negligible amount of scattering, the problem reduces to a rather simple 
situation. 

Let us suppose in a gross fashion that there are four contributions to the intensity I* (0, /x; r 3 ) 
in the direction /x at the top of the atmosphere (r = 0). One contribution will be from the atmos- 
phere above the cloud. The second will be from the cloud itself. The third will be from the at- 
mosphere beneath the cloud. The fourth will be from the ground. 
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Before formulating the problem a slight change in notation should be indicated. We shall in 
general indicate an intensity as 


X ( T a>M; r y -r p ) , 

where r a refers to the optical depth at which the intensity is considered, ^ refers to the cosine of 
the zenith angle of interest, and (r y - r ^ refers to the optical thickness of the layer in question. 

Thus I* (0, m; t 3 ) , the intensity of radiation in the direction ^ at the top of the atmosphere (r = 0) 
as indicated above, refers to the total optical thickness r 3 of the atmosphere- cloud composite*. 

We shall specify that the normal optical depths (r = 0 , r p r 2 , t 3 ) refer respectively to the top of 
the atmosphere, top of the cloud, bottom of the cloud, and the ground level. The azimuth- independent 
nature of the radiation field as implied by the lack of specification of the azimuthal angle & in our , 
notation should be evident in view of the fact that the problem has been stated in the context of * 
plane -parallel strata with no horizontal inhomogeneities. 

One point should be clearly emphasized; as the problem has been stated there is no scatter- 
ing of radiation except by the cloud itself. Thus aerosols, dust, and the ground are assumed to 
contribute no scattered energy to the radiation field. If this approximation is admitted in any 
particular case a considerable simplification results. We turn now to a formulation and solution 
of the problem of limb darkening in a cloudy atmosphere. Figure 4 illustrates the various inten- 
sity components which will arise in the course of our investigation. 

The nine assumptions specified in the Introduction will be maintained in the following discus- 
sion. In what follows we shall first want to formulate the expression for the intensity I c * {r l% ^ 
at the top of the cloud in the direction +pu Once 
the expression for I * (r 1? r 3 - has been 
obtained, we shall want to use this expression 
to formulate properly the expression for the in- 
) in the direction ^ at the top 


We proceed by considering first the con- 
tribution to the intensity I c * t 3 -rj from 

the atmosphere of optical thickness r x above the 
cloud. This component will be the intensity in- 
tegrated over all solid angles in the upward 
hemisphere which originates in the atmosphere 
above the cloud in the direction ( <f>') and is 

diffusely reflected by the cloud into the direction 
(m, <P)> Upon integrating over <// all the terms of 
the Fourier expansion in cos n ($>' -<f>) (cf. Equa- 
tions 152 and 157) drop out except for n = 0, 


tensity I* ( 0 , /x ; t 3 
of the atmosphere 



Figure 4 — The various intensity components arising 
through the interaction between an axially symmetric 
thermal radiation field and a model atmosphere- 
cloud-ground composite. The symbol ©refers to radia- 
tion originating in the relevant stratum. The tops of the 
upper atmosphere, cloud, lower atmosphere, and ground 
are situated respectively at r - 0, 7" lf r 2 , and t 3 . 
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leaving us with the expression 


j-f 1 


•>0 


T 1’“M 


r l ) w 


(181) 


where l A (r 1? -pi'; rj is the axially symmetric solution of the restricted atmosphere of optical 
thickness r x for the intensity at r x in the direction -/x\ 

Note that the intensity I A , -/x'; rj from the atmosphere above the cloud has no component 
due to the scattering of radiation in the atmosphere from a source (e.g., the cloud layer) outside 
this atmospheric layer. Of course radiation originating from the cloud, for example, is absorbed 
and re-emitted by the atmosphere above the cloud, and this will contribute to a change in the tem- 
perature profile in the atmosphere above the cloud. However, this temperature profile is presumed 
to be known (e.g., determined from radiosonde data) and thus i A (r 1? -yx'; r 2 ) known (for a given 
mixing ratio of absorbing constituents, etc.). 

The next contribution to the intensity I c *(t 1? /*; ~^i) to be considered originates from within 

the cloud itself. But this component of the intensity is simply 


1 2 ' I c ( T I- ^ T 2 _T l) ’ (l 82 ) 

the solution to the restricted problem (cf. Equation 117 and Appendix B). 

The third contribution originates in the atmosphere below the cloud, and will consist of two 
components. The first component is the intensity directly transmitted through the cloud and thus 
reduced by the factor e ^ . The second component is the intensity integrated over all solid 

angles in the lower hemisphere originating in the direction <£') and diffusely transmitted through 
the cloud into the direction (fi, 4>). 

These components are respectively given by 


I 


3 


'( t 2 -t i)A 


m ; r 


(183) 


and 


i 


4 



T 2 ~ T l) J A ( T 2 - M*; T 3 _T 2 ) <V • 


(184) 


Here the intensity I A (r 2 , /x; r 3 -r 2 ) is the restricted solution for the pertinent layer of atmosphere 
below the cloud. The transmission function T^°>(-ju. , -yx'; r 2 -rjJ is the azimuth- independent solution 
of the restricted problem for the "inverted" cloud. In other words, if we were to turn the cloud 
upside down and irradiate the cloud from what is now "above" in the direction -/x% the appropriate 
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transmission function becomes 


I c (r,, - M ; Tj-t,) 


4^ T c (0) (M. m': 



(185) 


The only reason for writing the transmission function as T c (0 )(-^, -/x' ; t 2 - t,) is to remind ourselves 
that the parity, or sense of direction, is reversed from what it is normally. 

The fourth and last contribution to the intensity 1 ^( 7 *!, r 3 is due to the isotropically 
radiating ground, here considered to be a blackbody radiating in accordance with Kirchhoff's law. 

This contribution also consists of two components. The first component is the intensity directly 
transmitted through the atmosphere below the cloud and thence through the cloud itself in the di- 
rection (m, 0). The second component is the intensity directly transmitted through the atmosphere 
below the cloud in the direction (m', 4 >')> and thence diffusely transmitted through the "inverted” 
cloud into the direction (m, 0 ). 

Denoting the intensity of the ground by the constant I G we have respectively for these two 
components: 
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“ ( T 3 ~ T 1 ) / ^ 
e / 


(186) 


and 
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Io e’ (T3 " Ta) / M * d M ' 


(187) 


where the integration in Equation 187 is once again implicitly performed over the lower hemis- 
phere, thus accounting for the notation of the transmission function. 

Upon adding up all the components Ij through i 6 indicated in Equations 181 through 187 we ob- 
tain for the outgoing intensity at the top of the cloud the expression 



The expression for the outgoing intensity I* (o, /z; t 3 ) at the top of the atmosphere may now 
be derived. This intensity may be thought to consist of two components. The first component is 
due to thermal emission in the atmosphere above the cloud alone and is nothing more than the 
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intensity corresponding to the restricted solution for such a layer of atmosphere. The second 
component is the intensity at the top of the cloud in the direction n as given by Equation 188 re- 
duced by the factor e Tl/ ^. These components are respectively 

I 7 = I A (0, W T,) (189) 


and 


T - X * / \ " T |A 

i b - : c ( T i * T i~ r i) e • (190) 

The exponential in Equation 190 is required to account for the reduction in intensity suffered by 
the radiation originating at the cloud top in traversing an optical path length (tj/m) through the 
atmosphere above the cloud. 

Upon collecting terms from Equations 188 to 190 we finally obtain for the outgoing intensity 
in the direction m at the top of the atmosphere the expression 


I* (°> MI t 3 ) = I* (0. M; Ti) + e Ti/M I c (t,, /x; 7 2 -t,) + e ^ I A (r 2 , /x; t 3 - t 2 ) 


+ T o e 3/ + ~^T 


(-/i, V: T a -Tj) I A ( 
+ s c (0) (m. m'; r 2 ~ r i ) i a( t i- V; 


( 7 2 1 ^ ; T 3 T 2j r ! G e 


- ( t 3 ~ t 2 )/ m ' 


(191) 


For a semi- infinite cloud Equation 191 reduces to 


I*(0,M) = I A (o. M ; Tj ) + e' Tl / w Ic ( TiiM ) 

e' T iA f 1 

+ 2^ s c <0) (m. m') 1*^, -m'; Tj) cV . (192) 

Jo 


Appendix D gives a convenient method for solving the integrals in Equations 191 and 192 by 
mechanical quadratures compatible with solutions to the restricted problems in the scheme of the 
n th approximation. In this case, we would seek solutions of the form 


I*(0,Ll.;t 3 ) * I* (0, M ; r 3 ) ( M = M . ; i = 1 , • • • , n) , (I 93 ) 

where are the intervals appropriate to n -point quadrature formula chosen for the purpose. 
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If we were to consider more than one scattering layer, such as two cloud decks for example, 
we would encounter a problem of somewhat more complexity than that of the first problem con- 
sidered in this section— the problem of describing the outgoing visible radiation field diffusely re- 
flected from a cloud-atmosphere composite. Since this would take us too far astray from our 
main theme, which is to discuss problems which are not excessively difficult to solve numerically 
on an electronic computer, we shall conclude our investigation at this point. 

The Albedo and Net Outgoing Flux 

Two quantities which follow immediately from the results of this section, where relevant, are 
the monochromatic albedo A* and net outgoing monochromatic flux. The net outgoing monochro- 
matic flux ttF follows from the relation 

F = 2 f fjL i*(o, m; T 3 )dM . (194) 

Jo 

The monochromatic albedo A* of visible radiation follows from the relation 

* - F 

A ^o F o ’ (195) 

where ttF 0 is the net monochromatic incident flux crossing a unit area normal to the direction of 
incident radiation. 

CONCLUDING REMARKS 

The analysis leading from the concept of individual photon-particle interactions to the final 
solutions for the angular distribution of outgoing radiation from the top of a cloudy atmosphere 
has required at various times the need for certain assumptions and approximations which are, to 
a greater or lesser extent, at variance with physical reality. At this point it would be well worth 
while to review briefly these various restrictions and, where possible, make a partial assessment 
of the limitations they will logically be expected to impose on the accuracy and validity of the 
solutions. 

In the first place, the polarization of the radiation field has been completely ignored throughout 
the formulation and solution of the equation of transfer. A proper account of polarization through 
the correct matrix formulation of the equation of transfer (References 3 and 9) is particularly dif- 
ficult and consequently has not been considered. However, most of the radiation scattered once is 
in general confined to a rather small range of scattering angles near zero degrees (Reference 10) 
and at © = 0° the degree of polarization must go to zero from considerations of symmetry. We 
would therefore expect the degree of polarization of the net radiation field always to be small, and 
the corresponding effect on the resultant intensity quite weak. For a clearer picture of how a 
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proper account of polarization affects the total intensity in the case of Rayleigh scattering, where 
forward scattering does not dominate and where the degree of polarization in general is consider- 
able, the reader is referred to Tables XV and XXIV and Figures 24 and 25 of Reference 3. 

The assumption that multiple scattering does not take place in a mass element dm which is 
large enough to contain a representative particle size distribution would appear to be quite sound 
for the wavelengths of radiation of interest. For a volume element of one cubic centimeter con- 
taining 300 particles, each having an efficiency factor for extinction Q E (i) = 3 and a radius r. = 20^ 
(microns) (References 6 and 11), it is seen from Equations 38 and 43 that, upon setting dz - l , 

dr % (300) (3) O) (20x 10“ 4 ) 2 cm' 1 dz % 0.01. (196) 

Since the numbers in Equation 196 are on the large side, it would be reasonable to assume that in 
general no more than a small fraction of a percent of any radiation scattered once in dm would be 
scattered again before leaving dm. 

In order to determine how valid the plane wave approximation of radiation incident on each 
particle is, we must obtain an estimate of how far an outgoing scattered wave must travel before 
becoming essentially a plane wave. It can be shown (Reference 5) that the magnitude of the radial 
components E r (s) and H r Cs > of the scattered wave are proportional to (A/r) 2 , where r is the distance 
from the scattering center and A is the wavelength of radiation, whereas the tangential components 
E t (s) and h t (s > of the scattered wave are proportional to (A/'r). If we impose the requirement that 
that magnitude of the radial components be about one percent of the magnitude of the tangential 
components before the outgoing scattered wave may be considered plane, we obtain the relation 

r > 100 A , (197) 

and from Equation 196 it may be inferred that more than 99% of all scattered radiation of wave- 
lengths less than A r 0. 1mm may be considered plane before suffering another scattering process. 

It should also be noted that, for particle densities of N 0 ^300 particles/cm 3 , the mean separation 
distance of particles is on the order of a millimeter, and this in turn implies that for wavelengths 
less than this order of magnitude all interactions of radiation with matter may be considered as 
"particle" phenomena; i.e., it is proper to consider clouds to be comprised of individual particles 
as opposed to a continuous fluid. No attempt will be made to discuss how deviations from spherical 
particles and/or particles containing impurities will affect the effective cross-sections or angular 
scattering patterns and the degree of polarization. 

The concept of local thermodynamic equilibrium in the pressure and temperature ranges of 
interest would appear to rest on an equally sound foundation for the atmosphere in general (Ref- 
erence 12). Since densities in the interior of particles are considerably greater than densities 
in the surrounding atmosphere, Kirchhoff T s law should be an even better approximation with 
regard to particle absorption and emission than for the surrounding atmosphere. 
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Treating the ground as an ideal blackbody in the infrared region of the spectrum is open to 
question. Variations of the emissivity of the ground from unity by several percent are quite pos- 
sible, although oceans would be expected to be quite black. 

The approximation of replacing an outside source with a point source removed to infinity was 
incorporated in the derivation of the equation of transfer which describes the diffuse radiation 
field. This approximation is also implicity extended to include any part of the radiation field wher- 
ever numerical solutions by quadratures ill the n th approximation are sought. The general validity 
of this method would appear to be born out from an inspection of results obtained by the author 
(unpublished) where the equation of transfer is solved in varying degrees of approximation for 
certain special problems and the results compared with one another. Little change was noticed in 
the solutions beyond the fifth approximation. By analogy with these results, it would also appear 
that any source extended over several degrees may be approximated rather accurately by a point 
source, so long as the extended source does not vary in intensity from point to point in too dis- 
continuous a manner. 

The replacement of the phase function for single scattering with an infinite series expansion 
of Legendre polynomials would be expected to yield quite rigorous results for any reasonably 
well-behaved phase function. The validity of approximating this function with a finite series of 
Legendre polynomials would, then, appear to depend only on the number of terms retained. In 
principle this is of no concern; in practice, however, the number of terms carried must be quite 
large, since the forward scattering nature of the phase function in general is quite extreme and 
the form of the function is not very amenable to reproduction by means of periodic functions. Un- 
published results by the author would indicate that the scattering and transmission functions of 
the previous section reflect the oscillatory nature of poorly fitted phase functions (low values of 
N), although mean curves fitted through the results appear fairly accurate within 2 % or so of the 
exact solution, at least for orders of N such that the magnitude of the oscillations in the approxi- 
mate phase function do not exceed 10% of the magnitude of the forward peak. For ratios of forward 
to backward scattering of 10 3 - 10 4 , the number of terms retained should correspond to N ~ 15-20 
or more. The solutions for the law of darkening do not tend to reflect these oscillations signifi- 
cantly, since multiple scattering and thermal emission play a much more dominant role, and the 
effect of single scattering is greatly reduced. 

A report by Neiburger (Reference 13) indicates that the particle size distribution in stratus 
clouds is fairly independent of optical depth except near the base of the cloud even though the 
liquid water content is quite dependent upon r. Whether this is true of clouds in general is doubt- 
ful, as is indicated from data obtained for a cumulus congestus cloud (Figure 48 of Reference 11). 
What is really required for thick clouds, however, is a knowledge of how constant the particle size 
distribution is near the boundaries of the cloud, since it is here that the contribution to the out- 
going diffuse radiation field is largest. More data, both observational and theoretical, are needed 
before a proper evaluation can be made of this problem. An inclusion of molecular effects from 
the pervading atmosphere compounds the difficulty. 
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The effect of scattering from a cubic centimeter of air at STP can be calculated from the ex- 
pression for Rayleigh scattering: 


N 


o * s 


8tt 3 ( n 2 - l) 2 
3 N 0 A« 


(198) 


Using the value of n = 1. 00029 , valid for 3000A, it can be verified that the optical thickness of 
a one centimeter column of atmosphere (assuming no absorption) is (Equation 43) 


dr 


( 10 ~ 


A 4 


N o X s dz * \ k / 


(199) 


where dz is set equal to one centimeter. If \^3000A, the optical thickness dr becomes about 


dr * l. 2 X 10 6 (200) 

Again, if we assume a particle density of N 0 = 100 particles/cm 3 , an efficiency factor for scatter- 
ing of Q s (i) = 2 , and a (rather low) average particle radius of r. = 5x 10“ 4 cm, it is found from 
Equation 199 that the optical thickness of a one centimeter column due to the particles alone is 


dr = (100) (2) (77) (5x 10~ 4 ) 2 % 1.57xl0 -4 , (201) 

or a factor of about 100 greater than the optical thickness resulting from molecular scattering. 

We conclude that scattering due to molecules contributes no more than one percent to the total 
scattered radiation at k = 3000 A, and this contribution becomes rapidly less as k increases. 

The final approximation to be considered is probably the most serious one— that of considering 
the whole problem in the context of plane-parallel layers containing no horizontal inhomogeneities. 
The really serious part of the difficulty with regard to real clouds is in connection with the in- 
herently "rough" surfaces of clouds being replaced in the theory by perfectly smooth and level 
surfaces. 

It is felt that no intuitive approach is going to be completely satisfactory. Intuitively one is 
led to believe that, in the first approximation at least, an inclusion of a rough surface in the prob- 
lem would lead to substantially the same results, since the roughness is only of a higher degree 
than that already incorporated into the problem implicitly through the concept of discrete parti- 
cles as opposed to the concept of a continuous fluid. Just as intuitively it might be argued that the 
outgoing radiation field will tend to become more isotropic so long as the roughness is shallow, 
because different orientations of the surfaces of each irregularity, and internal reflections among 
these irregularities, will tend to be random. If the roughness is great enough there will also be 
numerous shadow effects caused by each irregularity (although internal reflections among the dif- 
ferent individual irregularities will tend to reduce this effect), and the argument breaks down. 
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In this case one can even conceive of predominantly backward scattering, since the shadow effect 
is minimum in this direction. 

The reason the first intuitive approach is invalid is because the roughness of the cloud sur- 
face is not comparable to the "roughness” implicitly introduced into the equation of transfer through 
the concept of discrete scattering centers. Each mass element dm must now be considered large 
enough to contain a representative distribution of "irregularity” sizes, and, apart from considera- 
tions of the shapes, associated phase functions, and dependence upon optical depth of these ir- 
regularities, there is also clearly the impossibility of neglecting multiple scattering within dm it- 
self. On the other hand, if dm is restricted to a size such that only single scattering prevails, then 
clearly dm does not contain a representative particle size distribution, and is furthermore quite 
dependent upon real horizontal inhomogeneities. 

The other intuitive approach is invalid for the very reason that the whole problem was under- 
taken in the first place. In other words, the extreme anisotropic nature of the phase function for 
single scattering makes it impossible to discuss with any degree of confidence any tendency the 
radiation field might have toward isotropy at the surface of the cloud due to random orientations 
of the surfaces of the irregularities and internal reflections among them. 

If the task were not too prohibitive, a possible direction of investigation of these effects would 
be to consider the effect of different types of mass elements distributed in accordance with some 
distribution function, random or otherwise. For higher degrees of roughness it might be advan- 
tageous to pursue the problem from the point of view of considering each irregularity as an entity 
within itself, solve for the outgoing radiation field from each irregularity within the framework 
of an idealized geometry, and obtain the solution to the composite problem in a manner analogous 
to that adopted in the previous section. It is anticipated, however, that either approach is likely 
to be subject to a great deal more labor than the problem warrants. 

Three problems of quite general significance which have not been considered in this paper 
include: 

1. The dependence of the phase function for single scattering upon optical depth. 

2. The polarization of the radiation field. 

3. The solution to the equation of transfer within the framework of geometries other than 
plane- parallel. 

The first and third problems are of particular interest in the context of radiative transfer through 
aerosols, where the optical depth remains small for considerable linear distances. Furthermore, 
a correct account of polarization is of significant interest in connection with an interpretation of 
the characteristics of clouds in other planetary atmospheres, notably those of Jupiter and Venus. 
And finally, it is stressed that these problems should be attacked with a certain degree of sophis- 
tication. It is very difficult to deal intuitively with them, and the advent of high speed electronic 
computers has made the necessity for crude approximations less requisite than ever before. 
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Appendix A 


Normalization of the Phase Function for Single Scattering 


It is assumed that the phase function for single scattering, fP(cos0), can be expressed as a 
finite series expansion of Legendre polynomials in the form 


fP( cos 0) 



co t Pj (cos©) 


(Al) 


It is required to show that ?(cos©) expressed in this way is normalized to 2> 0 , the albedo for single 
scattering. This requirement is equivalent to the expression 

J *(cos0) ^ = 7L 0 , (A2) 

where the integration is performed over all solid angles. Thus 
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where spherical coordinates are chosen so that the polar angle is picked to be the scattering angle 
0, and the azimuthal angle <p is the angle of rotation about the axis of symmetry defined by 0 = 0. 

Upon reversing the order of integration and summation, and replacing cos ©by a, we obtain 




From the relations 
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we obtain 



and the required identity is established. 



Appendix B 




A Particular Integral to the Equation of 
Transfer Appropriate for an Emitting Cloud 


In context with the correct solution to the intensity distribution at the top of an emitting cloud, 
the equation of transfer appropriate to the restricted problem is 


dI(r, M ) _ 

M dr “ i( r ’ M) 


r+ 1 

- 4J pO> m') l(r, M 1 ) d/x' ~ (1 -w 0 )b(t) 


where p(^, ) has the form 
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P(M- M') = (m)P;. <m') 


i=o 


and b(t) has the form 


B(r) = £\ r r r 

r=0 


(Bl) 


(B2) 


(B3) 


A particular integral i[b(t)] satisfying Equation Bl will be required which, when added to the 
general solution to the homogeneous part of Equation Bl, gives the required solution for the angu- 
lar distribution of emitted radiation. 

Since we are concerned with solutions to the equation of transfer in the n th approximation, 
Equation Bl can be replaced with 2n equations of the form 
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From the manner in which it enters into the equation of transfer, and from the nature of the 
Planck function, we are led to assume that the required particular integral may be expressed in 
the form 


jBCr)] = ^ ^ £ C r,Ps-rK) ' 

r=0 s= r V 


(B5) 


where the C r S ’s are functions of a> l and b r only, and obey certain recursion relations to be derived 
below. 


From Equation B5 and the relation 


(21 +1)a*Pi(aO = (l +l)P i + 1 (M) 4 *Pi-i(M) 


Equation B4 becomes successively 
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the second term of the right-hand side of Equation B7 becomes successively 


= i|] a l p i(N5 Tr £ c r.sZ] a j p l (A 

= t T ' t p - ("0 ■ 

Therefore Equation B7 becomes, upon replacing r by (r + 1) on the left-hand side, 
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The left-hand side of Equation Bll is zero for r = s = N. Therefore 
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Since Equation B13 is valid for all r we require that 
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Equation B14 may be re-written as 
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for (o < r <N - 1). Note tha t the second term on the left-hand side of Equation B15 is zero for 
s = r - l; the summation is therefore started from s = r . 

Since Equation B15 must be valid for all /x. , we require that 


and 


(r +1) (N- r) 

2( N - r ) - 1 L r+i,N 


1 - 


2(N - r) + 1 


'r,N 


(0 < r <N- 1) 


(B16) 


(r + l)(N-r-l) 

2(N - r) - 3 r + i,N-i 


1 - 


^N-r-l 


2(N-r) -1 


(0 < r <N -2) 


(B17) 


Equations B16 and B17 were obtained by setting s = n and s = N - 1 respectively in Equation B15. 

From the identity P 0 = l, we further obtain from Equation B15, upon setting s = r, the 
relation 


3 


( r + 1) C 


r+1 ,r + 2 


(i 


for s = r(0£r <N-2). 

In the case s = r= N- lit follows from Equation B15 that 

^N“l,N-l ” ^N-l 


With the aid of Equations B12 and B16 through B19, Equation B15 reduces to 


N-2 N-2 

(r + l) 2(s-r)-l ( s-r ) P s-r(^i) + ( r + 1 ) 2(s-r)+3 (s~r + l) P,_ r (/i t ) 


(B18) 


(B19) 


N-2 



2(s ~r) + 1 




s= r+1 


Since Equation B20 must be valid for all M i} the general recursion relation 


1 


2(s - r) + 1 


(r + 1) 


(s ~ l + 1) 

2(s - r) + 3 '"r+1 ,s+ 2 


+ 


(S - r) 

2( s - r ) - 1 


C 


r + 1 , s 


(B20) 


(B21) 


for (0 £ r < N - 3; r + l<s<N-2)is obtained. 

In view of the fact that Equations B12, B16 through B19, and B21 are just sufficient to determine 
all the req uir ed values of C r s (o < r < N ; r < s <N ), we conclude that Equation B5 is a valid represen- 
tation for the required particular integral. 
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Two sets of values of C r s may be generated independently of each other, namely the set for 
which (N - s) is even and the set for which (N - s) is odd. Let us denote these sets by Equations 
B-I and B-n respectively. In what follows we shall assume N to be arbitrary. 


For Set B-I... 

From Equation B12 it is noted that 


*~N,N 


(B-I-l) 


From Equation B16 are generated successively 

N 

C N-1,N 7 1 ^ N C N,N 

l 1 " 3 “lj 
2(N-1) 


/ 1 ’Z^\ ^N-l, 

" I “2) 


L o,N 7 Y ^ l * N 

(2N-D ^nj 

To proceed, we obtain from Equation B18 


r - - . (N-D c 

Sf" 2 ,N -2 °N -2 T 3 ^ N-l.N 


From Equation B21 are generated successively 



(N 

"2) 

2 r , 

c 1 

3, N~2 

(*- 

1 - u 

3 *>1) 

5 C N-2,N 

^N-2 ,N“ 2 J 


(N 

"3) 

r3 c 

2 c 

‘4 , N~ 2 

(*- 

■5 “ 2 ) 

7 ^N“3,N 

3 U N-3,N~: 



1 

[>-i 

N - 2 

'O.N-2 

(, 

^N- 2 \ 

|_2N - 1 C i , n 

+ 2N- 5 


V 1 

2N -3/ 




Upon returning to Equation B18 we obtain 

bN - 4 + 3(1 ~£ 0 ) C N-3.N-2 


N -3 


(B-I-2) 


(B-I-3) 


(B-I-4) 


^N“4.N-4 


(B-I- 5) 
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Again, returning to Equation B21, we generate successively C N _ SiN _ 4 , C N _ 6 N _ 4 , • - • , C 0 N _ 4 . 
Thus Equations B18 and B21 are alternately employed as required until all the C r s ’ s for (N - s ) 
even are generated. 

The C r s ’s in Set (B-II), i.e., for (N - s) odd, are then obtained as follows. 

For Set B-II... 

From Equation B19 it is noted that 


^n-i,n-i kjv- 


(B-n-i) 


From Equation B17 are generated successively 


~N-2,N-1 


(N~ 1) 


/ 1 ^ \ S-1 ( N“1 

V ~ 3 "J 


^N-3,N“1 


2(N- 2) 




(B-II- 2) 


"O.N-1 


(N - 1 ) 


, N-l 

(2N-3) 11 - 


'l.N-l 


To proceed, we obtain from Equation B18 


c 

N~3,N- 

3 b N~3 

+ c 
3(1 -2t 0 ) L n-2.n-i ■ 

From Equation B21 are generated successively 

p — 

(N-3) 

r 2 r + „ 1 

N-4,N“3 1 

(■ - 1 a.) 

|_5 ^N-3, N-l Si-3, N"3j 1 

c = 

N-S.N-3 j 

(N-4) 

I" 3 

v 1 “ 5 S a) 

[7 C N~4,N-1 + 3 ^N“4,N“3 

p — _ 

1 

It" N - 2 N-3 

^O.N — 3 

' W N-3 \ 

2N -3 C i,n-i 1 2N - 7 C 

\ 

} 2N - 5 ) 


Upon returning to Equation B18 we obtain 


r 

V "N~5 , N*"5 

" V 5 + 

N-4 

3(1-S 0 ) C ™- 3 ' 


(B-II- 3) 


> (B-II- 4) 


j 


(B-n-5) 
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I 


Again, returning to Equation B21 we generate successively C N _ 6N _ 5 , c N - 7 , N - 5 , * * * , c o,n-s- 
Thus Equations B18 and B21 are alternately employed as required until all the C r s ’s for (N - s) 
odd are generated- The procedure is exactly the same regardless of whether N is even or odd. 

Since the C r s * s are uniquely determined, and there are no internal inconsistencies in their 
derivation, we conclude that (Equation B5) 


Ii[B(0] 


N N 

r-0 s =r 


is the required particular integral, where the C r s ’ s are determined in accordance with the pre- 
ceding discussion. 

In particular, for r = o, 


N 

I ; [B(0)] = 22 C 0 . s p s (m,) (i -- ± 1 , , ±n) , 

s=0 

and for r = r x , 

N N 

*i [®( T l)] = 2l T ‘ r Z] C r. S P s-r(Mi) (i = ±1, , ±n) 

r= 0 s= r 


(B22) 


(B23) 


Equations B22 and B23, along with the 2 n boundary conditions 1(0,-/^) = I (r, , ) = 0 

( i = 1 , • • • , n ) and the general solution to the homogeneous part of the equation of transfer, suf- 
fice to determine the appropriate law of darkening at the top of the cloud for the restricted problem. 


67 


I 



Appendix C 


The Flux Integral 


The net flux tt¥ v in the frequency range (i- , v + dv) is defined by 


7tF., 


n + 1 

Ml, M. 4 >) 

1 


(Cl) 


where the Integration is performed over all solid angles. If I v (r v , y, 0 ) can be written as 

N 

I„ (n> ( T „. M)cosn(0 o -0) 

the flux integral may be reduced to 


(C2) 


n = 0 


2^ mI v (0) (t„, m) 


d M - 


(C3) 


where the integration over the azimuth- dependent terms all go to zero. 


Diffuse Reflection and Transmission. 

If a single outside point source is considered to be the sole source of radiation, the relevant 
equation of transfer describing the intensity of radiation at any point in the cloud in any direction 
is, upon dropping the subscript v, 


dl(r, 0) 

^ &? 


I(t. m. <P) 



P(M» & 


y' , 0 ' ) I(r, yd , 0' ) dyu.' d0* 


1. 

4 


-r/MO 

e 


p(mi 4 >\ 


<*•<>) 


(C4) 
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The solution for l (0 > (r, h) for this equation of transfer, using the method of discrete ordinates, is 


I ^Me' kaT r_N_ 


1 ' 0> {r, Mi ) - ^ F °|/ 1 , k a + 


M. a e “ 


1 i k a 


CL~ 1 


IN 

1 ( _k -*) p *k) 


1-1 = 0 


e 


-t/mo r 


1 + Mi/Mo 




M = o 


(i = ±1, - , ±n) (C5) 


if Z Q / l and r 1 / oo. 77F 0 is the flux of incoming radiation crossing a unit surface area and normal 
to it. The net outgoing flux is found from Equations C3 and C5, after interchanging the order of 
summation, from 


a i p i(^i) 


F(t) = 2 2H 3i I<0) = kaT 7~!^ ^ 

i CL= 1 1 = 0 i 


i F ,r 


■"•'E 

£ = o 


l-^k a + 2 F »r,e ‘ < C6 > 

i = 0 i 


1 \ Y~» a i Mi P i(Mi) 


In deriving the recursion relation (cf. Equations 77-80) 


2k + 1 - 


W\+l 


k(\+ 1) “ A + 1 £\- 1 


(C7) 


use is made of the identity 


N 

(*) = €i ( x ) d a, t 


(x) (X = 0, ■■■ , N) 


1 = 0 


where D x j ( x ) is defined by 




a i P A (M j) P; (Mj) 
D A.I (x ) " 2 l+M iX 


(C8) 


(C9) 


Upon substituting \ = 1 in Equations C8 and C9, Equation C6 reduces to 


v iy 


f(t) = f ° \h m * e ' kaT ^ ^ +k ^ + £ m -* e+kaT ^ ( _k ^ + r » e ' rA ° 

La=l a - 1 


(CIO) 
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Upon setting \ = 0 in Equation C7, and remembering that cf_ 2 = 0, Equation CIO further reduces to 


/ — » V * 1 / +k T “k_T \ 

F(^) = " M a e ) 

l^a-l 


^0 ?0 e 


~ T A) 




y^. 




(C11) 


The net flux crossing the top of the cloud is found, upon setting r = 0, from 


F(0) = (1-S 0 )F 0 . 


_a= 1 


Mo r 0 


(C12) - 


ajnd crossing the bottom of the cloud, upon setting T ~ r v from 


F K) = (i-SoK- 


1 




La - 1 


“ Mo ^0 e 


r i Aol 


(C13) 


If the cloud is semi- infinite, all M_ a - 0 (a - 1, • • • , n), and Equation C12 becomes 


F(0 y - - {1-cZ q )F q < ir) + ^o>o 



(C14) 


i " to 
4 4* 


In the special case of conservative scattering (25 0 - l) it is clear that Equations C13 and C14 
are no longer valid relations. In this case, as can be verified by a direct substitution into the 
equation of transfer, the azimuth- independent term of the intensity in the case of a finitely thick 
cloud is given by (Equation 94) 


I (0) 1 T ’ Mi) 


4 F 0 ' 


“k T 

NI e a 


1 +Mj k a 


N 

^ l ( tk a ) P i (m i ) 


1—1 = 0 


M- a e a 


w 


L i = o 


e 


T /“o 


1 + Mi/M 0 


M M i ) + 3M ° { l 1 ” 3 “ l ) T + M i } M 0 +M n 


L j = o 


(C15) 


where M ia (a = 1 , • • • , n - l), M 0 , and M n are the 2n constants of integration. It is clear from Equa- 
tion Cll and the preceding discussion that in this case the flux integral must reduce to 


F(r) 



(C16) 
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Since 


a y 






a x^i " °: - 3 


(C17) 




Equation C16 reduces after some algebra to 

F(t) - (j. 0 F 0 M 0 


(C? 18) 


The flux integral F (r) is seen to be independent of optical depth, as it must be for conservative 
scattering. 




. s 


If the cloud is semi- inf inite and = l, Equation C18 must reduce to 


F(r) = m 0 F 0 ; 


(0 19) 


i.e., the "albedo" of the cloud must be unity. The expression for I (0) (t, m £ ) in this case is seen 
by inspection to be 


I(0) ( T ’ Mi) 


NL e 




4 0 1/ , 1 + M i k a 

a = 1 


€ l ( + k a) (Mi ) 


I— t= 0 


e 


r Ao 


1 + Mi/Mo 




L 1 = o 


3m 0 


(l~3 X,) r +Mi +M n J 


(C20) 


where M a (a = 1, • • • , n - l) and M n are the n constants of integz^ation to be determined from the n 
boundary conditions 


y X v 




ic°) (0, -Mi) = 0 (i = 1 , •• * , n) 


(C21) 


The Law of Darkening. 

The equation of transfer appropriate to the case of axial symmetry (no outside sources) where 
thermal emission is taken into account is (Equation Bl) 


M dr 


:(w) 


P„0*. M') l(r v . M ')«^' - fl-S 0 )B„(T) 


(C22) 
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Upon dropping the subscript v, the solutions for I (r, fx) for this equation of transfer, using the 
method of discrete ordinates, are 


v’x e " k “ T r 

c ( T ’^i) = / rrjnr a 2_ s *^( +k .) p »(^) 

( 1 Lj = o 


y-M. a e +k “ T 

k « 


- w 


Li = o 


N 

r = 0 


n 

Z>.. p ~h) 


(i = ±1, ■ • • , ±n) . (C23) 


The relation for the net flux is therefore (cf. Equations C5, C6, and Cll) 


F(t) -• 4(l-S 0 )y iT" ( M -a e+k “ T ” M a e" k “ T ) + 2 


r=0 


N 

E c '- p - w 


(C24) 


Upon rearranging the order of summation, the second term on the right-hand side of Equation C24 
becomes successively (cf. Equation C17) 


2 / a. 


, 'b' W ' y P s-r( M 0 r - 2/ T S / / > ^ r. s P l(Mj) p s -r(^,) 


7/ C, rtl r p . (C25) 


The final expression also results from the identity 


where 





28 


n, m 

2m + 1 


f 0 (n / m)l 

\l (n = m)J 


(C26) 


(C27) 


Upon collecting the results from Equations C24 and C25, the expression for the net flux 
becomes 


F (r) = 4 


(^o) (n 


v , N~1 

')'iL 


(C28) 
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The expressions for the net flux at the top and bottom of the cloud are respectively 


F(0) 


n 

(l-"o) J]iT( M -a- M a) + JC 0i1 


(C29) 


and 


r ( T 0 


= 4 


n 


+k r. 


k t \ 1 

) + 3 


N-l 

L 

r = 0 


" r, r + 1 1 


If the cloud is semi- infinitely thick, Equation C28 reduces to 


(C30) 


n N-l 



a. - 1 r = 0 


and Equation C29 reduces to 


(C31) 


F(O') 


n 


- 4 


(1 -o) 



a = 1 


\ 

3 



(C32) 
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Appendix D 


A Note on the Solutions of Certain 
Integrals and Integral Equations 


Consider the integral 


f 


gO> m’) hO'> x") 


G(m. m") 


(Dl) 


If g ( m, ) and h (//. fi" ) correspond to so-called n th approximation solutions to the restricted 

problem, the integral in Equation Dl may be solved by replacing the integral with the sum 


G (/l; , fJy ) 


n 



j=l 



(D2) 


for discrete values of m, and /x k . The divisions are the positive roots of the Legendre polynomial 
P 2 n (m) (0 </x< 1 ), and the weights a j may be determined from the formula (Equation 61) 


d/x 


_i_r 

■ wft I, 


/X - /X 


d/x . 


Now consider the integral equation 


>(m. M 0 ) “ f (M» M 0 ) + I gO> At') h(/x', /x") S(/x\ ^ 0 ) d^u" d/x' 

Jo Jo 


(D3) 


(D4) 


where f, g, and h are known functions and S is to be solved for. 
Let (Equation D2) 


GO, /x") 


r 1 

gO» AOh(/x\ /x") d i/x' 

Jo 


% a. g( M , Mj ) h( Mj , K ") 

j = » 


(D5) 
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Then Equation D4 becomes 


*(m> M 0 ) " f ( m> M 0 ) + G(m. m") S(m"> M 0 ) d M" . 

Jo 


(D6) 


or, 


s(m, M 0 ) * f (M. M 0 ) + 


n ” 

22 a > g ^’ h ( M J' m ") 


S(m". M 0 ) d M" 


n , 

f (M,Mo)^f k a. 

k= 1 ^ = 1 


>(m> Mj) h( Mj , Mk) 


s(/^ k , m 0 )/ 


(D7) 


Not let m take on the discrete values 


M “ Mj (i - 1, • • • , n) 


(D8) 


Equation D7 now becomes 


s (Ml- m 0 ) - f (Mi. M 0 ) + 



■ n 

J2 a i g tv M i) h (*V 


S(M k . Mo) 


(D9) 


We have from Equation D9 a set of n equations in n unknowns of the form 


s (Mj. Mo) - t (mj > M 0 ) + ^"l, j ,i S (mi i Mo) + ‘ ‘ + Cjj n S(^ n , m 0 ) 


S(m 2 > Mo) f(M 2 ’Mo) + ^- 2 ,j,i S(mj , M 0 ) + ' ’ ' + ^2,j, n S(m„> M 0 ) 


(DIO) 


S (m„.M 0 ) - f (M n - Mo) + C nj i S^J, /I 0 ) + • • • + C n >j>n s(m„, Mo) 


where 


C. . . a, 

i , ) , k k 


- n 

Y2 a > g ( Mi> ^) h b’ 


Li = i 


( i = 1, • • • , n ; k = 1, , n) . (Dll) 
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We thus may obtain from Equations DIO and Dll the n values S^, m 0 ) which correspond to the 
function s(/x, fi Q ) through the expression 

s(/x, ,u 0 ) ^ /X 0 ) (fj. = /z.; i = 1, •• • , n) . (D12) 

The reader should be reminded that the solutions of Equations D2 and D9 depend implicitly in 
the present context on the solutions to the restricted problems which in turn require the evalua- 
tion at times of certain functions which are undefined for the divisions 1 = * * * > n )* 

More explicitly, it is required to evaluate the quantity 


where 


V = h<-»(m 0 )h(">(-m 0 ) 


(D13) 


H< m) (x) 


Mi 


£l ( X + Mj) 

j=l 

ftfl+ka"*) 

n = 1 


(D14) 


in the n ,h approximation, and then evaluate the functions 

•y m 

= 1 - VMo ‘ (D15) 

We note that Equation D15 is undefined for y. i - . Equation D15 may, however, be evaluated 

quite easily by writing it in the form 


3 :(m) (^ i .M o) - -m 0 (m 0 + mJ Mo) -Mo) 

where H (m > (pt., x) is defined by the expression 




ft (* * 


n ( x+ Mj) 


i»i+l 


n ( 1 + k a mx ) 


a'l 


(D16) 


(D17) 


Equation D16 is defined for jjl 0 in the range (of> 0 < l) including the values m 0 = Ab (* “ 1 > * ' ’ > n )> 
but excluding the values (m 0 = l/k a m )(cf. Equation D17). It is never necessary, however, to include the 
latter values of jj. 0 in any theoretical discussions of outgoing intensity distributions since other 
more convenient values of are just as useful. 
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It should be noted that, for low values of n, Equation D3 will in general not be the most ad- 
vantageous one to use in determining the appropriate weights needed in evaluating Equations D2 
and Dll, and the roots (mj , M k ) of P 2n (m) are correspondingly the incorrect discrete divisions 
of evaluation- The integrations in Equations D1 and D4 are over the interval (0</i<l), and Equa- 
tions D2 and Dll evaluate these integrals over the interval (“1 1 m 1 1). In effect this is accomplished 
by implicitly setting the relevant functions in the integrands identically equal to zero in the interval 
(-1 <ii <o); however this has the effect of causing the integrands to be discontinuous at (/x = 0) and 
therefore impossible to approximate accurately with a polynomial in ^ of low degree (small number 
of terms). 

A better method (Reference 8) of evaluating Equations D1 and D4 is through the use of the 
expressions 


l 

a i - T a i (D18) 

and 

= T( 1 + ^i ) (i = ±1, • •• , ±n) - , (D19) 

yielding the 2n values a. and ~jl i which are to replace the n values a. and/x A (i - l, * • * , n) in 
Equations D2 and Dll. This artifice effectively maps the interval (-1 < m < 1) onto the interval 
(OlMll). 

In order to reduce the number of terms involved in the relevant equations, it should suffice 
to solve now for the roots (/x.) of P 2n (/x) for low values of n. it will also be required to develop 
some kind of interpolation scheme for obtaining the values of functions at the divisions and 
which were previously determined numerically at other divisions for higher orders of n. The ef- 
ficacy of the method will be most manifest where it is desired to restrict the number of computa- 
tions to a minimum while still obtaining accurate results. For details of this method the reader 
is directed to the reference previously cited. 
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Appendix E 


List of Symbols 


The list of symbols is divided into the following seven categories: Coordinates; General 
Symbols; General Functions; Subscripts; Superscripts; Energy Functions; and Intensity, Scatter- 
ing, and Transmission Functions. 

For the most part this list is intended to be a guide and is in no sense complete. In many 
cases only the form of the symbol is given with the understanding that the explicit meaning of each 
symbol is to be obtained from the text; this requirement is expedient in view of the fact that to 
give a complete list, a great deal of repetition would be involved, and this would be both a waste of 
space and wearing upon the reader. The division of the list into various categories was performed 
with the aim of isolating certain symbols, such as scattering and transmission functions, which 
the reader will have recourse to most frequently. 

Parentheses enclosing a number in the right column refer to an equation in which the symbol 
is defined or illustrated. 


Coordinates : 

b, b‘ 

M» m' . m" 

Mi , Mj, M k . Mi 

0 , <&" 

$ o 

r ’ T v 


Zenith angle 

Cosine of the zenith angle treated as a continuous variable 
Cosine of the zenith angle treated as a discrete number 
Cosine of the zenith angle of an outside point source 
Azimuthal angle 

77 radians removed from the azimuthal angle of an outside point 
source 

Monochromatic optical depth 
Geometric height 


General Symbols: 

A, A' Area 
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A* 


a i’ a j ■ a k> a i 


C i.i.k 


D i 

E 

F 

F o 

g 

H 


I 

5 

k, k a , k a m 
M±a* M ± m a 
m 
N 


n 

a. 

n 


"l 

No 

P, (1 > (j = S,A) 


Albedo (195) 

Gaussian weights (61) 

Coefficients of the terms of a power series expansion for the 
Planck function (B3) 

Coefficients of scattering functions involved in a system of linear 
algebraic equations (Dll) 

Coefficients of the terms in a series expansion for a particular 
integral (86-91) 

Constant of integration (54) 

Particle size distribution (16) 

Energy (122); electric field vector 

Related to net flux (Cl) 

Related to incident flux (46) 

Probability of extinction by particles (122) 

Magnetic field vector 

Monochromatic specific intensity (1) 

Rate of emitted radiation (127) 

Roots of characteristic equations (79, 104) 

Constants of integration (76, 98) 

Mass 

Number of terms in series expansion (63, 69) 

Degree of approximation (70); real part of refractive index 
Complex refractive index 

Number of particles per unit volume per unit radius range cen- 
tered about r. (12) 
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Total number of particles per unit volume (16) 

Probability of scattering (S) or absorption (A) of one photon by a 
particle of radius r ; in dm 


Qj (i) (j = E, S, A) 


r 

S 

T 

t 

V 

x 


a 


P 


0 


K 

k 


V 


P 


a 


T a (a / v) 

Yj (j = E, S, A) 


CO 


CO 


0 


CO 


m 

l 


Efficiency factor for extinction (E), scattering (s), or absorption 
(A) of a particle of radius r. 

Radius; distance 

Surface (Figure 2) 

Temperature 

Time 

Volume 

Argument of function; in general indicates that the function is 
one of several arguments 

Solid angle; variable of integration 

General argument of function 

Angle through which radiation is singly scattered (Figure 1) 

Imaginary part of complex index of refraction 

Wavelength; running number of summation 

Frequency of radiation 

Density 

Area 

Normal optical thickness 

Effective cross-section for extinction (E), scattering (S), or 
absorption (A) of mass element dm (38) 

Solid angle 

Effective albedo for single scattering (39) 

Coefficients of terms in associated Legendre polynomial series 
expansion of the phase function for single scattering 


General Functions : 
B V (T)1 
B (r) J 


Dl, A ( x )| 
D U(*>J 


Planck function (29, 69) 
Characteristic functions (78, 103) 
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f(M) 

f (M, Mo) 1 

f (n) K = M. M 0 )J 

G(/x. M')\ 

G( Mi . M k )J 

g(M. M') ") 

g (n) ( T t ; m. m') J 

g(<£) 

H(x) j 
H< m > (X)J 

K (m> (m ; , x) 

Km'. M ") 
h i <n> (^- *o)l 

h 2 (n) (m. M")J 
h(r) 

J( T * 

0< n) (/x, Mo ) 1 
0C") (t [ ; M. Mo)J 
Pj (x) 

Pf (X) 

p(/i, m'. 0') 


f?( cos ©) 



8(a-/3) 
(x) | 

€?(*) J 


Function associated with H functions (D15) 
Undetermined function (47) 

Inhomogeneous terms of integral equations (168, D4) 
Integrals (Dl, D2) 

Factors in integral equations (169, D4) 

Undetermined function (47) 

H functions (84, 108) 

Modified H function (D17) 

Factor in integral equation (D4) 

Integrals (163, 164) 

Undetermined function (50) 

Source function 

Integrals (176, 180) 

Legendre polynomial 

Associated Legendre polynomial 

Phase function for single scattering (4, 6) 

Phase function for single scattering (3) 

Products of H functions (83, 107) 

Kronecker delta function (110, B9) 

Dirac delta function (47) 

Xi functions (80, 105^ 


Subscripts :* 

A Absorption; atmosphere 

C Cloud 

♦ All numerical subscripts except (0) refer in general to explicit values or functions. 
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D 

E 

G 

i 


j 

k 

l 

M 

m 

n 

P 

r 

S 

s 

T 

z 

a 

P 

k 

v 

0 


Diffuse 

Extinction; emission 
General; ground 

Discrete value or function; explicit particle size; running number 
of summation 

Discrete value; running number of summation 

Discrete value; running number of summation 

Degree of function; discrete value; running number of summation 

Molecule 

Discrete value 

Degree of function; discrete value 
Particle 

Degree of function; running number of summation; radial 
Scattered 

Degree of function; running number of summation 

Transmitted; tangential 

Height 

Discrete value; running number of summation 
Discrete value 

Degree of function; running number of summation 
Frequency 

General subscript— meaning to be derived from text 


Superscripts: 

i Explicit particle size; discrete value 

k 1 Associated with or explicit order of function or argument; order 

m > of azimuth- independent terms in the Fourier series expansion 

n J of a function 
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0 


Zero - order 


Energy fun c tions : 
E v (z, /z, 0) 


SE^ ( z, (i, 0) 
d[X(z, M . *)] 


Energy at a level (height) z which crosses an element of area in uni 
time and is restricted to the frequency range (v, v + dv) and to 
an element of solid angle containing the direction (as 0). The 
level, frequency, and directional dependence in general may or 
may not be indicated 

That fraction of E v (z, ji 9 0) which is restricted to the time in- 
terval (t, t + dt ) 

That fraction of SE V (z, as 0) which has arisen through an inter- 
action with matter in the mass element dm 




That fraction of E v (z, /i, 0) which is extinguished by the matter 
in dm 


d[AE v ( z t fi ' , 0' )] 


That fraction of AE V (z, 0) which is scattered only into the 

solid angle d co' containing the direction (m', <P * ) 


Intensity , Scattering , and Transmission Functions : 

Because of the great number of intensity, scattering, and transmission functions found in the 
text, only the different forms will be indicated in the list of symbols. With the aid of the ac- 
companying definitions this reduced list of symbols should suffice to allow a ready interpretation 
of any such function from the form of the function, even though the specific notation will vary. 


Intensity functions of the form. . . 


1(T, /X, 0) \ 

I v (z> 0)/ 

I A (T, fl, 0) 


Intensity of monochromatic radiation at respectively an optical 
depth r or a height z in the direction (as 0) 

Same as I (r, as 0) above referred to a restricted solution for a 
stratum of gaseous atmosphere 


I c (r, A X. 0) 


Same as I (t, as 0) above referred to a restricted solution for a 
cloud 


I A * ( r > M. 0) 
I c * (T, /x, 0) 


Same as I A (r, as 0) above referred to a composite solution; e.g., 
a component due to a neighboring cloud is added to i A (r ? as 0) 
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Same as I c (r, as 0) above referred to a composite solution; e.g., 
a component due to a neighboring stratum of gas is added to 

l c ( r > th 0) 



I* (0, M , <*>) 
I ( T a> m; r y - 

i a ( t <* 1 m ; 

J c ( T - MS T y 
I c ( T «* m; t, 

I* (o, M ; Tj) 

! c (w) 

I*(0, M) 

I< n) (T, At) 

I c 

li [ B (^)] 

It H 

Io <0) (-. Mi) 
d[s ± I. (z, /X, 


Resultant intensity at the top of the atmosphere (r = o ) for the 
total composite problem 

Tg) Intensity of monochromatic radiation at a normal optical depth 

r a in the direction m referred to a stratum of normal optical 
thickness ( T r ~ T y3 ) 

- r^) Same as I (r a , fi; r y - above referred to a restricted solution 

for a stratum of gaseous atmosphere of normal optical thickness 

( T r" T /3) 

Same as I [r a9 r y - above referred to a restricted solution 
for a stratum of cloud of normal optical thickness (r y 

,-Tp) Same as l c (r a , \x\ r T y ~ r ^) above referred to a composite solu- 

tion; e.g., a component due to a neighboring stratum of gas is 
added to I c (r a , m; * r p) 

Same as l*(0, m, 0) above in the direction m instead of (n, 0) re- 
ferred to an atmosphere- cloud composite of total normal optical 
thickness r 3 

Same as l c (r a , yx ; r y - r ^ above for a cloud of infinite optical 
thickness 

Same as I* (0, yx; r 3 ) above for r 3 infinite 

The azimuth- independent factor of the n th term in the Fourier 
series expansion of I (r, yx, 0) above 

Intensity of monochromatic thermal radiation isotropically 
emitted by the ground (186) 

A particular integral of the equation of transfer appropriate to 
a restricted solution to the problem of limb darkening (B5) 

A particular integral of the equation of transfer appropriate to 
a restricted solution to the problem of diffuse reflection and 
transmission (82) 

The general solution to the azimuth- independent equation of 
transfer 

0)] ( j = S, A, E) A loss or gain [respectively (-) or (+)] to the intensity of the 
radiation field at z in the direction (m, 0) due to scattering (S), 
absorption (A), or emission (E) processes with particles in the 
radius range (r., r. + drj contained in the mass element dm 
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S ± I. (z, /x, 0 ) (j - s. A, E) 


Scattering functions of the form, 

S(r ^ /lx, 0 ; /x' , 0 ') 

S A (tj; /x, < P \ m' - 0') 

S c (/x, 0; fji‘ , &' ) 

s a ( T i ; _ m> 0; - m . ) 

s c* ( T i ; ^ 1 ‘O 

s c (0 > (m, m') 
s c (0) (m. m' ; t 2 -Tj) 

S A (n > (At. /*') | 

S A (n) ( T l- M’ m')J 

S c (n ) (M, m') 

s c *(m>m # ) 1 

S c ( T i’ M. M') J 
S A (n) (-/x, V) 

S * (t x ; m , 0; M 0 . 0 O ) 

S* (n) (M> Mo) | 

S * (n) (r^ M. A^o) J 


The same as d [8 ± i. ( z , /x, 0)] above integrated over all particle 
sizes 


Scattering function (cf. Equation 129) appropriate to a stratum 
of normal optical thickness r x for radiation scattered from the 
direction (-/x', 0') into the direction (0, 0) 

Same as S(r x ; 0; /x', 0'] referred to the restricted solution 
appropriate to a stratum of gaseous atmosphere 

Same as S(r 1 ; /a, 0; m', 0 ') above referred to the restricted solu- 
tion appropriate to a semi- infinitely thick cloud 

Same as S a (tj ; /x,0; /x', 0 ') but in the reverse sense, i.e., the 
reverse side of the stratum 

Same as S c (/x, 0 ; /x # , 0' ) above referred to the composite prob- 
lem; e.g., a contribution is added to S c (/x, 0 ; /x',0' ) due to an 
overlying atmospheric stratum of normal optical thickness r x 

The azimuth-independent analog to S c (,u, 0 ; m', 0') above 

Same as S c CO) (m, m') above for a cloud of optical thickness (r 2 

Azimuth- independent factor of the n th term in the Fourier series 
expansion of S A (r 1 ; 0 ; /x' , 0 ') above 

Azimuth- independent factor of the n th term in the Fourier series 
expansion of S c (/x, 0; /x', 0 ' ) above 

Azimuth- independent factor of the n th term in the Fourier series 
expansion of s * (tj ; /x, 0 ; /x' , 0') above 

Azimuth- independent factor of the n th term in the Fourier series 
expansion of S A -/x , 0; -/x' ? 0 ') above 

Scattering function appropriate to the solution of the total com- 
posite problem at the top of the atmosphere 

Azimuth- independent factor of the n th term in the Fourier series 
expansion of S*(r 1 ; /x, 0; /x 0 , 0 O ) above 
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Transmission functions of the form . . . 


T(tjJ m» <f>\ 0') 

t a ( t i ; Mi 0 ; 0 ') 

t a ( ^ i ; ~m', 0') 


Transmission function (cf. Equation 130) appropriate to a 
stratum of normal optical thickness r x for radiation originally 
in the direction (— /x*, 0' ) which is transmitted diffusely through 
this stratum into the direction (— /x, 0) 

Same as T^; m, 0; m', 0') referred to the restricted solution 
appropriate to a stratum of gaseous atmosphere 

Same as T A (r x ; m, 0; m # , 0') but in the reverse sense, i.e., the 
reverse side of the stratum 


T a ( T 1 » Ml & m' . 0') 


T A (n > (Mi M') 


T A (n) (-M. -M') 

t c (0) (m. m'; T 2 ” T i) 


T (°) I- 

C 1 




T 2" T l) 


Same as m, 0; m', 0 # ) above referred to a composite solu- 

tion; e.g., a component due to a neighboring cloud is added to 
T A ( T 1 5 m,0;m'> 0') 

Azimuth- independent factor of the n th term in the Fourier series 
expansion of T A ; m, 0; m', <£') 

Azimuth- independent factor of the n th term in the Fourier series 
expansion of T A (rj 0; V, 0') 

The azimuth- independent transmission function appropriate to 
the restricted solution for a stratum of cloud of optical thickness 

( T 2- T l) 

Same as t c <°>(m, m'; t 2 _r i) above except in the reverse sense 
through the cloud 


NASA-LangLey, 1965 G- 595 
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"The aero?iautical and space activities of the United States shall he 
conducted so as to cofitribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere aiid space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof." 

— National Aeronautics and Space Act of 1958 
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